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Morphological  instability  refers  to  the  tendency  towards  spatial 
pattern  formation  on  the  liquid-solid  interface  when  a  dilute  binary 
mixture  is  solidified  or  fused.  The  importance  of  this  phenomenon  is  in 
the  growth  of  metal  alloy  and  semiconductor  crystals  from  their  melts, 
where  it  influences  the  solute  or  dopant  concentration  resulting  in 
nonuniform  physical  and  electrical  properties. 

Previous  formulations  of  morphological  instability  have  involved 
several  simplifying  assumptions  which  restricted  it  to  the  study  of  a 
region  immediately  surrounding  the  interface.  The  models  have  limited 
validity  and  they  require  separate  treatments  for  different  situations 
like  freezing  and  melting.  In  this  study  a  new  uniform  approach  is 
presented  which  considers  the  entire  melt  and  crystal  domain  and  is 
applicable  to  all  situations.  Earlier  formulations  are  shown  to  be 
approximations  of  this  and  exchange  of  stabilities  is  proven 
asymptotically. 


vx 


The  model  is  then  used  with  a  weakly  nonlinear  technique  to  predict 
the  shape  of  the  bifurcation  diagram  for  various  cell  patterns.  The 
subcritical  nature  of  morphological  instability  is  shown  and  regions  of 
its  prevalence  are  determined  over  the  entire  domain  of  experimental 
parameters.  This  was  used  to  compare  with  experimental  results  and  to 
determine  optimal  crystal  growth  regions. 

A  comprehensive  comparison  of  morphological  instability  with  con- 
vective  instability  was  undertaken  and  this  phenomenon  was  shown  to 
resemble  Marangoni  convection  in  its  mathematical  and  physical  fea- 
tures. This  was  done  in  order  to  introduce  some  of  the  multitudinous 
mathematical  techniques  employed  in  convective  instabilities  into  mor- 
phological instability  and  specifically  was  used  here  to  complete  the 
eigenspace  of  the  linearized  problem. 

Two  imperfections  which  reside  in  the  domain,  heat  loss  at  the 
container  wall  and  advection  in  the  melt,  were  considered  and  shown  to 
be  bifurcation  breaking  imperfections.  Solutions  to  the  problem  were 
obtained  in  both  cases  by  matched  asymptotic  expansions  and  based  on 
these  results  a  practical  way  of  minimizing  the  effect  of  these  imper- 
fections was  suggested. 
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CHAPTER  1 
DESCRIPTION  OF  THE  PROBLEM 


Morphological  instability  refers  to  the  process  of  spatial  pattern 
formation  at  the  liquid-solid  interface  when  a  binary  mixture  is 
solidified  or  fused.  This  is  a  problem  of  hydrodynamic  instability  and 
like  all  other  problems  of  this  nature,  for  this  phenomenon  too  there  is 
an  onset  point  where  the  initially  planar  interface  first  begins  to 
deform  and  forms  cellular  patterns.  These  grow  into  deeper  finger-like 
shapes  and  eventually  forming  side  branches  and  tree-like  dendritic 
structures. 

The  importance  of  this  phenomenon  is  in  the  growth  of  metal  alloy 
and  semiconductor  crystals  from  their  melts.  Morphological  instability 
affects  not  just  crystal  shape  but  also  the  solute  or  dopant 
concentration,  resulting  in  nonuniform  physical  and  electrical 
properties  in  the  crystal.  This  is  especially  perfidious  in 
applications  where  superfine  crystals  with  very  consistent  properties 
are  required.  Recently  there  have  been  some  indications  that  crystal 
quality  can  be  significantly  improved  by  growing  them  in  low  gravity  as 
this  reduces  other  problems  associated  with  crystal  growth  like  natural 
convection,  but  unfortunately  not  morphological  instability.  Growing 
the  crystal  at  very  high  temperature  gradients  or  at  very  low  growth 
rates  avoids  morphological  instability  but  crystals  grown  at  high 
temperature  gradients  are  of  poor  quality  due  to  thermal  stresses  and 
the  very  low  growth  rates  makes  the  process  very  expensive.  Hence  this 


becomes  a  problem  not  only  of  avoiding  crystal  surface  deformations  but 
one  of  optimization  of  the  process  as  well. 

There  are  several  methods  for  growing  crystals  from  the  melt, 
distinguished  by  the  hydrodynamics  of  growth.  The  three  basic  ones  are 
Bridgman,  Czochralski  and  float  zone  and  most  techniques  are  variations 
of  these,  like  horizontal  and  vertical  Bridgman.  A  typical  Bridgman 
experiment  is  shown  in  Fig.  1-1.  The  material  is  usually  in  a  quartz 
ampoule  and  is  melted  and  then  recrystallized  in  the  Bridgman  furnace. 
The  upper  part  of  the  ampoule  is  maintained  at  a  higher  temperature  than 
the  lower  and  solidification  proceeds  upwards  and  the  ampoule  is  pulled 
downwards  at  the  same  velocity  v.  The  top  of  the  melt  is  protected  by  a 
liquid  encapsulant  like  B2CU.  At  the  end  of  the  process  the  ampoule  is 
broken  to  retrieve  the  crystal. 

In  the  float  zone  technique  shown  in  Fig.  1-2,  the  ampoule  or  the 
material  itself  is  pulled  through  a  circular  furnace  and  the  melting  and 
recrystallization  proceeds  simultaneously.  As  this  technique  can  be 
done  even  without  a  container  it  avoids  the  problem  of  impurities  from 
the  ampoule  entering  the  crystal,  but  its  more  difficult  to  maintain  a 
uniform  temperature  gradient.  Figure  1-3  shows  the  Czochralski  method 
where  the  crystal  is  rotated  and  pulled  from  a  melt  pool.  As  the 
emphasis  here  is  on  the  liquid-solid  interface,  the  modelling  of  the 
crystal  growth  process  will  be  kept  as  general  as  possible  but  will 
resemble  the  Bridgman  technique  the  most. 

The  temperature  and  concentration  profiles  during  typical  crystal 
growth  conditions  are  shown  in  Fig.  1-4.  The  temperature  profiles  in 
the  liquid  and  solid  are  virtually  straight  lines  and  solute 
concentration  in  the  solid  is  virtually  constant.    But  the  solute 
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concentration  in  the  liquid  C}  changes  sharply  near  the  interface 
because  of  solute  rejection  on  solidification.  This  in  turn  has  an 
effect  on  the  freezing  point  depression  in  the  liquid  as  shown  in  the 
figure.  Figure  1-5  shows  the  same  profiles  but  in  a  situation  where  the 
freezing  point  in  the  liquid  TM  now  exceeds  the  actual  liquid 
temperature  (this  change  can  be  brought  about  by  either  reducing  the 
liquid  temperature  gradient  or  by  making  the  change  in  C.  near  the 
interface  even  sharper  by  increasing  the  growth  velocity) .  This  is 
referred  to  as  constitutional  supercooling  and  the  system  responds  to 
this  unstable  situation  by  interface  deformation.  Countering  this  is 
the  interfacial  tension  which  always  acts  to  minimize  the  surface  area 
which  in  this  case  is  the  planar  surface.  When  this  balance  is  upset, 
or  in  other  words  when  the  onset  conditions  are  exceeded,  the  interface 
loses  planarity  and  forms  cellular  patterns.  Figure  1-6  shows  the 
profiles  for  the  fusion  case,  where  we  could  have  TM  in  the  solid  being 
less  than  the  actual  solid  temperature  and  once  again  interfacial 
deformation  is  the  system's  response,  balanced  by  interfacial  tension. 
The  only  difference  here  is  that  since  solid  diffusivities  tend  to  be 
several  orders  of  magnitude  lower  than  liquid  diffusivities,  the  solute 
concentration  profile  in  the  solid  near  the  interface  will  vary  even 
more  sharply  resulting  in  lower  onset  conditions  for  morphological 
instability. 

In  the  paper  of  Trivedi  and  Somboonsuk  (1984)  there  is  a  series  of 
photographs  from  an  experiment  (see  Fig.  1  in  their  paper)  where 
succinonitrile/acetone  crystals  were  grown.  The  first  photograph  shows 
the  liquid-solid  interface  just  after  onset  and  has  a  discernible 
cellular  pattern.   Later  ones  show  the  cells  becoming  deeper,  forming 
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Figure  1-4.   Concentration  and  temperature  profiles  during 
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Figure  1-6.   Concentration  and  temperature  profiles  during  fusion 


fingers  and  ending  up  as  dendrites.  The  arrows  in  the  first  photograph 
mark  the  initial  perturbations  that  eventually  become  dendrites.  In 
this  thesis  we  will  concentrate  only  on  the  region  near  the  onset  point, 
shown  by  the  first  two  photographs,  though  dendrites  will  be  mentioned 
in  discussions. 

In  the  papers  by  Morris  and  Winegard  (1969)  and  Tiller  and  Rutter 
(1956)  we  see  another  aspect  of  morphological  instability,  the  variety 
of  cellular  patterns,  finger-like  shapes,  hexagonal  cells  and  variations 
of  these.  Other  possible  shapes  are  cylindrical  rolls  and  rectangular 
cells  but  by  far  the  commonly  encountered  pattern  is  the  hexagonal 
one.  The  choice  of  the  cell  pattern  is  extremely  important  and  factors 
determining  this  choice  will  be  discussed  later.  Their  figures  also 
show  that,  unlike  other  forms  of  hydrodynamic  instability,  the  number  of 
cells  on  a  single  crystal  is  in  the  hundreds. 

Though  it  is  customary  to  model  morphological  instability  in  terms 
of  the  temperature  and  concentration  profiles,  in  reality  this 
phenomenon  seldom  exists  in  isolation;  it  is  usually  coupled  with  fluid 
flow  in  the  melt.  There  are  two  kinds  of  flows  that  occur.  The  first 
is  buoyancy  driven  solutal  convection  which  is  caused  by  the  sharp 
solute  concentration  gradients  in  the  melt.  The  other  is  the  result  of 
density  change  during  solidification.  When  solidification  occurs  there 
is  a  constant  rate  of  volume  decrease  which  causes  the  melt  to  move  in 
to  fill  the  vacated  space.  This  motion  is  referred  to  as  advection  and 
the  rigid  side  walls  of  the  ampoule  will  cause  closed  streamlined  flow 
in  the  melt  as  a  result  (see  Fig.  6-2).  In  addition  there  will  be  flows 
in  the  melt  in  Czochralski  growth  due  to  rotation  and  other  kinds  of 
flows  in  special  growth  techniques. 


Apart  from  these,  several  other  parameters  affect  this  phenomenon, 
the  most  important  of  which  are  due  to  the  fact  that  most  crystals  are 
faceted;  that  is,  they  have  a  crystal  lattice  structure.  Hence  whether 
the  lattice  axis  is  aligned  or  not  with  the  growth  direction  is 
extremely  important  as  can  be  seen  from  the  experiments  of  Heslot  and 
Libchaber  (1985).  Other  important  considerations  are  grain  boundaries, 
wetting  of  the  ampoule  wall  and  the  presence  of  impurities.  Also  in 
rapid  solidification,  the  system  will  not  be  at  thermodynamic 
equilibrium  and  kinetic  undercooling  of  the  melt  becomes  significant. 


CHAPTER  2 
PREVIOUS  WORK  ON  MORPHOLOGICAL  INSTABILITY 


2.1 .   Early  Work 

The  first  successful  attempt  at  explaining  morphological 
instability  qualitatively  was  by  Rutter  and  Chalmers  (1953).  They 
coined  the  word  "constitutional  supercooling"  to  describe  the  existence 
of  unstable  melt  regions  near  the  interface  where  the  freezing 
temperature  can  be  higher  than  the  liquid  temperature  itself  and 
correctly  identified  this  as  the  cause  of  interface  deformation. 

Tiller,  Rutter,  Jackson  and  Chalmers  (1953)  quantified 
constitutional  supercooling  and  for  instability  came  up  with  the 
condition 


mGic   >  "  Gl  (2.1-1) 


where  m  is  the  absolute  value  of  the  liquidus  slope,  G?  the  liquid 
temperature  gradient  at  the  interface  and  G.  the  solute  concentration 
gradient  in  the  liquid  at  the  interface.  The  negative  sign  is  caused 
by  G.  and  G.  being  in  opposite  directions  (see  Fig.  1-5).  As  can  be 
seen  these  simple  thermodynamic  explanations  do  not  take  into  account 
the  stabilizing  effect  of  interfacial  tension.  To  do  so  would  require 
casting  the  problem  as  one  of  hydrodynamic  instability  and  obtaining  the 
onset  conditions  from  a  linear  stability  analysis.  This  is  exactly  what 


Mullins  and  Sekerka  (1963,  1964)  did  when  they  considered  the  problem 
with  temperature  and  concentration  equations  in  the  liquid  and  solid  and 
boundary  conditions  at  the  interface.  Their  criterion  for  instability 
to  an  infinitesimal  disturbance  was 


mG„  >  -  GT  +  a2TMX/L.  (2.1-2) 

c      T      M   h 


where  GT  is  the  weighted  temperature  gradient 


gt  "  (k*Gz  +  ksGa)/(ka  +  ks}       (2-1"3) 


Here  kg  and  k.  are  the  solid  and  liquid  thermal  conductivities.  G  is  a 
modified  liquid  concentration  gradient  given  by 


Gc  =  GlQUl  "   v/2V/(a*  "  (1/2  "  k)v/V         (2.1-4) 


2     2    2  1/2 
a£  =  (a"  +  y/HDp1^  (2.1-5) 


where  a  is  the  wavenumber  of  the  disturbance,  k  the  solute  distribution 
coefficient,  TM  the  melting  temperature  of  the  pure  solvent,  Lh  the 
latent  heat  of  fusion  and  \  the  interfacial  tension.  This  analysis  laid 
the  foundation  for  all  further  work  in  morphological  instability. 
Following  them  Woodruff  (1968)  did  the  linear  stability  analysis  along 
the  same  lines  for  the  melting  problem  and  came  up  with  the  same 
criterion  as  (2.1-2)  but  with 


G  =  G  (na  -  v/2D„)/(nka  +  a„  +  (1  -  k)v/2Dn)        (2.1-6) 

C      SC    S         X,        S      X,  X, 
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2     2    2  1/? 
a  =  (a*  +  vV4D:),/<:  (2.1-7) 

s  s 


where  Ggc  is  the  solid  concentration  gradient  at  the  interface  and  n  is 


D  /D.,  the  ratio  of  solute  diffusivities. 

S   jC 


2.2.   Later  Research 

The  next  major  contribution  to  the  problem  was  made  by  Wollkind  and 
Segel  (1970)  who  proved  "exchange  of  stabilities"  for  this  problem  for 
most  parameter  ranges.  Proving  exchange  of  stabilities  is  equivalent  to 
showing  the  existence  of  the  onset  of  steady  state  nonplanar 
solutions.  They  also  considered  the  weakly  nonlinear  regime  after  the 
onset  of  instability  and  using  the  method  of  Stuart  (1960)  and  Watson 
(1960)  analyzed  the  problem  for  the  case  of  two-dimensional  rolls, 
showing  the  existence  of  subcritical  bifurcation  at  most  growth 
velocities. 

The  method  of  Stuart  and  Watson  is  essentially  the  theory  of  Landau 
(see  Drazin  and  Reid  (1981))  and  follows  the  dominant  mode  of 
instability  into  the  weakly  nonlinear  regime.  In  this  form  it  is 
applicable  only  to  disturbances  of  one  cellular  pattern  at  a  time,  but 
Segel  and  Stuart  (1962)  extended  this  theory  to  the  prediction  of  the 
preferred  pattern  by  considering  the  interaction  of  two  specified  modes 
of  disturbance.  Depending  on  the  way  these  two  modes  were  combined  it 
was  possible  to  obtain  two  dimensional  rolls  or  hexagonal  cells  and  they 
showed  that  the  experimental  parameters  would  determine  the  stability  of 
these  patterns.  Sriranganathan,  Wollkind  and  Oulton  (1983)  adopted  this 
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method  for  morphological  instability  and  gave  parameter  ranges  where 
each  type  of  cell  was  stable.  The  limitation  of  the  method  is  that  it 
considers  hexagonal  and  two  dimensional  roll  patterns  but  not 
rectangular  cells  or  cylindrical  rolls  and  ignores  the  effect  of 
container  shape  and  size  which  have  been  shown  to  be  important  in  wave 
pattern  selection  (see  Koschmeider  (1967)). 

Ungar  and  Brown  (1984a)  considered  the  highly  nonlinear  problem  and 
after  making  several  simplifications  obtained  solutions  using  the  finite 
element  method.  Finite  elements  can  handle  highly  nonlinear  problems 
and  give  very  accurate  numerical  solutions  but  are  extremely  time 
consuming.  Solving  the  full  morphological  problem  is  a  very  expensive 
proposition  by  this  method  and  hence  Ungar  and  Brown  simplified  the 
problem  by  ignoring  the  latent  heat  and  solid  diffusivity  and  assuming 
that  thermal  conductivities  in  liquid  and  solid  were  equal.  This 
allowed  them  to  reduce  the  problem  to  a  "one-sided  model"  consisting  of 
variables  in  the  liquid  region  only,  considerably  simplifying  the 
algebra  and  saving  computer  time.  Such  a  model  will  have  only  limited 
validity  in  highly  nonlinear  regions  and  this  was  borne  out  when  Ungar, 
Bennett  and  Brown  (1985)  solved  the  complete  problem.  But  their  most 
extensive  calculations  were  done  only  for  the  one  sided  model  and  hence 
this  is  of  chief  interest.  These  were  done  only  for  the  case  of  two 
dimensional  roll  disturbances  and  here  they  showed  that  contrary  to  that 
reported  by  Wollkind  and  Segel ,  there  were  multiple  regions  of 
supercritical  and  subcritical  bifurcation.  More  importantly  they  showed 
that  at  large  deformations  of  the  interface  secondary  bifurcations 
occurred.  Ungar  and  Brown  (1985)  also  modelled  the  formation  of  deep 
cells  in  an  attempt  to  follow  the  transition  to  dendritic  growth. 
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Nonlinear  finite  difference  calculations  were  done  in  a  more 
limited  way  by  McFadden  and  Coriell  (1984)  for  the  two  dimensional 
case.  Later  McFadden,  Boisvert  and  Sekerka  (1987)  extended  the 
calculations  for  the  three  dimensional  patterns  of  hexagons  and  cross 
rolls.  In  both  cases  the  enormous  expenses  involved  restricted 
calculations  to  a  few  parameter  values. 

2.3.  Inclusion  of  Other  Effects 

While  these  workers  were  investigating  the  basic  problem  others 
were  busy  trying  to  incorporate  various  influences.  The  most  important 
concern  was  the  effect  of  fluid  flow.  Delves  (1968,  1971  and  1974)  in 
attempting  to  approximate  the  influence  of  advection  and  stirring  in  the 
melt,  studied  the  influence  of  plane  Couette  flow  on  the  problem.  He 
showed  that  two  dimensional  roll  disturbances  in  the  flow  direction  were 
stabilized  but  there  was  no  effect  on  disturbances  perpendicular  to  the 
flow.  Coriell,  McFadden,  Boisvert  and  Sekerka  (1984)  modelled  Couette 
flow  more  systematically  and  came  to  the  same  conclusion.  Recently 
McFadden,  Coriell  and  Alexander  (1988)  examined  the  effect  of  plane 
stagnation  flow  on  two  dimensional  disturbances  perpendicular  to  the 
flow  and  here  too  the  flow  as  found  to  be  stabilizing. 

In  another  very  important  development  Coriell,  Cordes,  Boettinger 
and  Sekerka  (1980)  studied  morphological  instability  with  solutal 
convection.  They  showed  that  the  two  instabilities  were  essentially 
decoupled  with  the  melt  being  unstable  to  convective  disturbances  of 
long  wavelengths  and  the  interface  unstable  to  nonplanar  disturbances  of 
small  wavelengths.   Also  at  low  growth  rates  the  dominant  instability 
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was  convective  and  the  interface  was  not  easily  disturbed.  At  high 
growth  rates  the  roles  were  reversed  and  at  an  intermediate  velocity  the 
two  instabilities  became  comparable.  It  was  only  at  this  rate  the  two 
instabilities  interacted  and  the  result  was  the  prevalence  of 
oscillatory  instabilities.  Their  conclusion  was  that,  except  at  this 
particular  growth  rate,  it  is  usually  sufficient  to  study  only  the 
dominant  instability  near  its  onset. 

Following  Coriell  et  al .  several  workers  have  looked  at  special 
aspects  of  these  two  instabilities  and  their  work  has  been  reviewed  by 
Glicksman,  Coriell  and  McFadden  (1986).  They  all  confirmed  or  refined 
the  work  of  Coriell  et  al.  but  all  the  main  conclusions  mentioned  above 
still  hold. 

Several  other  influences  apart  from  fluid  flow  have  been 
incorporated  into  the  model  but  only  a  few  relevant  ones  will  be 
considered  here.  Coriell  and  Sekerka  (1972,  1973)  tried  to  include  the 
effect  of  grain  boundaries  on  morphological  instability  by  assuming  that 
its  only  effect  was  to  shift  the  onset  conditions.  They  failed  to 
observe  that  in  the  presence  of  grain  boundaries  there  could  be  no 
planar  solutions  to  the  problem  and  that  the  interface  will  be  nonplanar 
at  all  times.  Ungar  and  Brown  (1984b)  obtained  the  solutions  to  this 
problem  by  matched  asymptotic  expansions  for  small  grain  angles  and 
using  finite  elements  for  solutions  of  large  grain  angles. 

In  rapid  solidification  kinetic  undercooling  of  the  melt  is 
significant  and  Seidensticker  (1967)  included  this  and  showed  that  it 
caused  a  shift  in  the  onset  conditions.  The  significance  of  this  was 
shown  by  Hardy  and  Coriell  (1968,  1969  and  1970)  when  they  observed 
morphological  instability  in  the  growth  of  ice  crystals.  Constitutional 
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supercooling  was  not  a  factor  here  and  it  was  shown  that  kinetic 
undercooling  was  the  primary  cause.  This  dual  cause  for  morphological 
instability  is  somewhat  analogous  to  the  situation  in  natural  convection 
where  we  find  that  the  variations  of  density  and  surface  tension  with 
temperature  can  both  cause  convective  instability. 

2.H.     Experiments  in  Morphological  Instability 

The  early  work  on  modelling  morphological  instability  was  prompted 
by  experimental  observations  but  beyond  that  very  few  quantitative 
experiments  have  been  done  near  the  onset  conditions.  This  is  an 
unfortunate  state  of  affairs  and  experimental  verifications  of 
theoretical  predictions  are  badly  needed  if  further  concrete  progress  on 
the  theoretical  front  is  to  be  made.  The  work  of  Morris  and  Winegard 
(1969),  Trivedi  and  Somboonsuk  (1984)  and  of  Heslot  and  Libchaber  (1985) 
have  already  been  mentioned.  Recently  de  Cheveigne,  Guthman  and  Lebrun 
(1985,  1986)  have  attempted  to  verify  the  weakly  nonlinear  and  strongly 
nonlinear  theoretical  predictions  and  one  hopes  that  more  work  along 
these  lines  will  follow. 

De  Cheveigne  et  al .  performed  their  experiments  on 
succinonitrile/acetone  and  CBr4/Br2  systems.  (These  organic  mixtures 
are  much  easier  to  work  with  than  metal  alloys  as  they  are  generally 
nonfaceted,  transparent  and  require  small  temperature  gradients  and 
hence  they  have  been  very  popular  with  experimentalists.)  They  found 
that  the  cell  pattern  formed  and  its  dimensions  were  strongly  dependent 
on  geometry  of  the  container.  More  importantly  when  they  ran  the 
experiments  for  two  dimensional  roll  patterns,  they  found  only 
subcritical  instability. 
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2.5.  Limitations  of  Existing  Models  and  Unaddressed  Issues 

In  Chapter  1  the  cause  of  constitutional  supercooling  was  explained 
as  being  due  to  the  sharp  solute  concentration  gradient  in  the  liquid 
near  the  interface,  while  elsewhere  in  the  liquid  and  the  solid  the 
solute  concentration  was  practically  a  constant.  It  would  seem  then 
that  the  only  region  of  interest  is  the  interface  and  a  liquid  "boundary 
layer"  adjacent  to  it.  This  has  prompted  all  previous  workers  in 
morphological  instability  to  consider  D  /v  as  the  characteristic  length 
of  the  problem  and  to  ignore  solid  diffusion.  A  typical  value 
of  D^/v  is  100  microns  and  this  means  that  the  far  ends  of  the  melt  and 
crystal  are  infinitely  far  away  and  the  domain  of  the  problem  is 
effectively  confined  to  the  liquid  boundary  layer  mentioned  above.  For 
the  melting  problem  a  characteristic  length  of  D_/v  is  used  and  the 
domain  becomes  an  even  smaller  boundary  layer  in  the  solid. 

These  assumptions  considerably  simplify  the  algebra  involved  and 
hence  their  popularity.  But  they  constrain  the  validity  of  the  model  in 
several  ways.  The  most  obvious  one  is  that  they  necessitate  the  melting 
and  solidification  problems  to  be  studied  separately,  even  though  they 
only  differ  in  the  direction  of  the  growth  velocity.  Besides  this 
assumption  fails  for  very  small  growth  velocities,  as  it  introduces  a 
singularity  at  v  =  0.  Later  we  will  show  that  neglecting  solid 
diffusion  also  introduces  a  singularity  and  makes  the  model  fail  in  the 
nonlinear  regime. 

Finally,  any  effect  which  resides  in  the  entire  domain,  not  merely 
the  boundary  layer,  cannot  easily  be  incorporated  into  the  model,  which 
is  why  all  influences  on  morphological  instability  studied  so  far  are 
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either  boundary  layer  effects  (e.g.,  solutal  convection)  or  interfacial 
effects  (e.g.,  grain  boundaries  and  kinetic  undercooling).  Phenomena 
that  span  the  entire  domain,  like  advection  in  the  melt  or  imperfect 
insultion  of  the  ampoule  walls,  have  been  either  inadequately  treated  or 
ignored  completely.  Hence  there  is  a  need  for  a  model  that  includes  the 
entire  liquid  and  solid  domains  which  would  be  applicable  for  all  growth 
velocities.  This  model  should  also  dispense  with  the  separate 
treatments  accorded  so  far  to  the  solidification  and  fusion  problems 
with  one  uniform  formulation. 

In  Section  2.2  it  was  mentioned  in  connection  with  the  work  of 
Wollkind  and  Segel  (1970)  and  of  Ungar  and  Brown  (1984a),  that  this 
problem  oscillates  between  subcritical  and  supercritical  instabilities 
for  the  case  of  two  dimensional  roll  disturbances.  They  did  not, 
however,  compute  the  ranges  of  each  type  of  instability  for  the 
experimental  parameters  involved.  This  is  necessary  in  light  of  the 
experiments  of  de  Cheveigne  et  al.  (1986)  who  observed  no  supercritical 
instability.  Also  the  extension  of  these  predictions  to  three 
dimensional  disturbances  like  hexagonal  and  rectangular  cells  is  yet  to 
be  done. 

It  would  not  be  an  exaggeration  to  state  that  the  inspiration  for 
all  the  theoretical  work  done  so  far  in  morphological  instability  has 
come  from  Rayleigh-Benard  convective  instability.  A  comparison  between 
the  two  problems  would  be  invaluable  as  a  source  of  continued 
inspiration  and  as  a  way  to  draw  conclusions  and  conjectures  about 
morphological  instability  from  the  vast  published  literature  on 
Rayleigh-Benard  convection.  Hurle  (1985)  has  attempted  this  but  his 
work  can  only  be  regarded  as  perfunctory  and  there  exists  a  need  for  a 
more  rigorous  treatment  of  the  issue. 


CHAPTER  3 
A  UNIFORM  FORMULATION 


3.1.  The  Formulation 


Since  we  are  not  making  the  assumption  that  the  liquid  and  solid 
are  very  deep,  the  problem  has  to  be  formulated  very  carefully, 
especially  with  regard  to  the  outer  boundaries,  if  we  are  to  avoid  an 
intractable  moving  boundary  problem. 

A  typical  crystal  growth  set  up  is  shown  in  Fig.  3-1 .  The  ampoule 
is  heated  by  the  heating  coils  surrounding  it  and  they  keep  the  melt 
region  at  a  temperature  T1  and  the  crystal  at  Tg.  The  temperatures  T-, 
and  T2  are  maintained  constant  by  means  of  thermocouples  located  at  z  = 
s  and  2  «  -%.  The  ampoule  is  pulled  towards  the  cooler  end  at  the  same 
velocity  V  at  which  the  crystal  grows,  thus  keeping  the  interface 
stationary.  The  region  near  the  interface  is  protected  by  an  insulating 
shield  and  it  is  this  region  that  becomes  the  domain  in  our  model. 

So  in  this  model  the  outer  liquid  and  solid  boundaries  become  fixed 
at  z— I  and  z=s  respectively  and  the  solid  will  be  moving  with  a  bulk 
velocity  v  and  the  liquid  with  a  bulk  velocity  v/Y,  where  Y  is  the  ratio 
of  densities  Pg/p.*  In  this  section  we  will  assume  that  Y=1  and 
consider  the  effects  of  Y  not  being  unity  in  Chapter  6  as  this  would 
cause  advection  and  fundamentally  alter  the  basic  problem.  Also  we  will 
assume  that  the  melt  concentration  at  the  outer  liquid  boundary  is  a 
constant  C,  . 
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Figure   3-1.      Experimental  set  up 
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The  domain  equations  in  the  liquid  melt  are 


dTl  dh  „2  <3.1-1 ) 


3t  +  v  II  =  V  h 


act    ac^  (3.1-2) 


8t  +  V  TI  =  V  Cl 


In  the  solid  region  the  equations  are 


3T      9T  (3.1-3) 

—2.  +  v   s  .  0  v^T 
dt      9z    s   s 

9C      9C  (3.1-4) 

-sl  +  v  -jS.  =  D  V^C 

9t      9z    s   s 

where  T  and  C  are  the  temperature  and  solute  concentration,  D  and  a  the 
diffusivity  and  thermal  diffusivity,  with  the  subscripts  I     and  s 
referring  to  the  liquid  and  solid. 
The  boundary  conditions  are 


T4  "  V  C4  "  C1  at  z  =  ~l  (3.1-5) 


Ts=T2  at  z=s  (3.1-6) 


At  the  liquid-solid  interface  we  will  use  x,   to  denote  the  departure 
from  planarity  and  write  the  boundary  conditions  at  z=z 


h  =  Ts  =  TM-mC£-TMr-"  (3-1~7) 

h 


kAVT*  *  "  -  ksVTs  *  n  =  Lh(v  +  DT}  (3-1"8) 
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kC^   =  Cg  (3.1-9) 


VC*    •»-(▼*  dT)C£   =  VCs    •  a  -   (v  ♦  §i)Cfl  (3.1-10) 


where  TM  is  the  melting  point  of  the  pure  solvent,  A  the  interfacial 
tension,  L^  the  latent  heat,  m  the  absolute  value  of  the  liquidus 
slope,  k.  and  k  the  liquid  and  solid  thermal  conductivities,  k  the 

•C       s 

distribution  coefficient,  n  the  normal  at  the  interface  directed  into 
the  solid  and  H  is  the  curvature  of  the  liquid-solid  interface  and  in 
Cartesian  co-ordinates  is  given  by 


n   L  .  2  V1    v3x;  ;     3x  3y  3x3y   .  2  ll    v3x;  n 
3x  3y 


[1  +  (|f)2  +  (|f)2]"3/2  (3.1-11) 


In    cylindrical    co-ordinates,     if    we    assume    circular    symmetry,     it 
becomes 


H-   [     £\  (1    ♦   (|)2)   Ah    *   (g)2]_3/2  (3.1-12) 

3r 


It  is  assumed  that  the  side  walls  are  sufficiently  far  apart  and 
well  insulated  to  enable  us  to  impose  periodic  boundary  conditions  in 
that  direction. 

To  convert  these  equations  into  the  dimensionless  form  we  will  use 

the  liquid  depth  I   as  the  characteristic  length  and  the  diffusive  time 

2 
D?/v  as  the  characteristic  time. 


21 


The  temperature  will  be  made  dimensionless  by  f  =  (T-T  )/G  I   with 

M   T 


GT  "-  (k*G*  +  ksGs)/(ks+  V  (3'1~13) 

where  G^  and  Gg  are  the  temperature  gradients  in  the  liquid  and  solid 

in  the  quiescent  planar  state. 

Similiarly  the  concentration  will  be  made  dimensionless  by 

C  =  C/G  I   with 
c 

Gc  "  (V*c  +  DsGsc)/(DJl  +  Ds}  (3-1"l24) 

where   G   and  G   are  the  concentration  gradients  in  the  liquid  and 
solid  in  the  quiescent  planar  state. 

So  in  dimensionless  form,  if  we  neglect  the  Lewis  numbers  D./ct  in 
the  temperature  equations,  we  have 


V2T£  =  0  (3.1-15) 


3C    _  3C  (3.1-16) 

-A   +  v  _*  .  v2 

3t      3z      * 

V2Ts  =  0  (3.1-17) 


3C    _  3C  (3.1-18) 

— ■  +  v  — §■  =  nV^C 
3t      3z       s 

where  v  =      and  n  =  D  /D0 

The  boundary  conditions  are 


at  z=  -1 , 
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h    -    (T1  -  TM)/GTl  "  T1 


ct  -  va0i  -  c, 


(3.1-19) 


at  z=s/£=s, 


Ts  =  (T2  "  W  "  T2 


(3.1-20) 


at  the  interface  z=£ 


\-\-  'SeCZ  ~A  « 


(3.1-21) 


7T  •  n  -  6VT  •  n  =  L(v  +  2i) 

*         S  Dt 


(3.1-22) 


kCfl  =  C 
I         s 


(3.1-23) 


VC„  •  n  -  (v  +  ^)c0  =  nVC  •  n  -  (v  +  ^)C 

1  Dt  %  S  Dt   S 


(3.1-24) 


where  —  is  the  total  derivative,  6  the  ratio  of  thermal  conductivities 

ks/k£  and  Se'  ^  and  L  are  the  Sekerkai  capillary  and  latent  heat  numbers 
respectively. 


Se  = 


A  = 


mG 
c 

M 


L  = 


w 

LhD£ 
k,GT£ 


(3.1-25) 


(3.1-26) 


(3.1-27) 
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At  this  point  a  discussion  of  the  choice  of  a  critical  parameter 
becomes  imperative.  The  experimentally  variable  parameters  for  this 
problem  are  G^,  C  and  v  or  their  equivalents  in  this  formulation 
GT'  Gc  and  v"  Most  Previous  workers  have  adopted  G.  or  C  as  their 
critical  parameter  but  Hurle  (1985)  has  proposed  Se,  by  analogy  with  the 
Rayleigh-Marangoni  problem.  Recently  de  Cheveigne  et  al .  (1986)  have 
advocated  the  use  of  v  from  an  experimentalist's  perspective.  Although 
this  is  a  valid  choice  the  reason  no  one  else  has  used  it  so  far  is 
probably  because  v  occurs  in  the  domain  equations  and  will  give  rise  to 
an  infinite  number  of  eigenvalues  in  the  linearized  problem.  We  second 
Hurle's  suggestion  and  choose  Se  as  this  seems  to  be  the  naturally 
occurring  coupling  factor  between  temperature  and  concentration  and  the 
fundamental  cause  for  morphological  instability.  Besides  it  includes 
both  G  and  G  .  But  the  suggestion  of  de  Cheveigne  et  al .  still 
remains  a  valid  one. 

3.2.  The  Linear  Stability  Problem 

From  now  onwards  the  "  -  "  for  dimensionless  variables  will  be 
dropped.  The  steady  state  planar  solution  to  this  problem  occurs  when 


C  (x,y)  or  r  (r)  =  0  (3.2-1) 


The  solution  is  then 


Tlc(z)  =  -  (T1  +  SeCac(o))z  -  SeC^o)      (3.2-2) 
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Tsc(z)  =  (T2  +  SeCiQ(o))z/s   -  SeC^Co)  (3.2-3) 


C   (z)  =  C  r   d-k)exp(vz)      i  r  (l-k)exp(-v)       i       ,    ., 
°)lcu;   c1Lk-exp(-v(l+s/n))   1J/Lk-exp(-v(1+s/n))   1J      (3<   } 

c  fzl  -  r  r  (l-k)exp(vz/n)    ,i/r   (l-k)exp(-v)      i      ,     . 
CscU)  "  C1Lkexp(v(1+s/Ti))-1   1J/[k-exp(-v(1+s/n))   1J      (3"2  5) 

To  write  the  equations  of  the  linear  stability  problem  we  will 
impose  an  infinitesimal  disturbance  on  the  steady  state  solution. 


h  '  ho   +  h  (3'2-6) 


with  similar  expansions  for  T  ,  C.,  C  and  £. 

9    m  S 

Considering  the  linear  stability  problem,  in  order  to  separate 
variables  we  will  assume  a  horizontal  cell  pattern.  This  pattern  for 
two  dimensional  rolls  is  given  by 


i    \        r-        2irnx        ...  ,  2irn  ,  _    _   „, 

<j>   (x)   =  Cos  — : — ,   with  wave  number     a     =  — —  (3.2-7) 

n  L  n        L 


*.  i-i  /        \  2imx        „        2imy  2irn  ,,   .   0, 

for  cross  rolls     $   (x,y)   =  Cos  — —  +  Cos  — ^-  ,      a     =  —r—  (3.2-8) 

n  u  u  n         L 


for  rectangular  cells 


4   (x,y)   =  Cos  Zl™  cos  2EHC,   a     =  2im(1/L2   +   1/L2,)  V2  (3.2-9) 

n  Li .  i_*  _  n  i  £- 


for  hexagonal  cells 
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♦„(».»>  -  2Cos  ffS  Cos  2E*  .  Cos  ^.  an  -  <g      (3.2-,0, 


for  cylindrical  rolls  $  (r)  =  J  (a  r/R)  (3.2-11) 


where  an  are  the  zeros  of  J1 ,  and  JQ  and  J.  are  the  zeroth  order  and 
first  order  Bessel's  functions  of  the  first  kind.  We  can  now  write 


T£(x,y,z,t)  =  Tz1(z)  4>](x,y)eat  (3.2-12) 


with  corresponding  forms  for  the  other  variables,  where  a  is  the 
eigenvalue  of  the  linear  stability  problem. 

The  linear  stability  problem  becomes  in  the  domain 


(D2  -  a2)Tz1=  0  (3.2-13) 


oC£1=  (D2  -  vD  -  a2)Cn  (3.2-14) 


(D2  -  a2)Tg1=  0  (3.2-15) 


iS,  =  (D2-^-a2)C  ,  (3-2"l5) 

n  si        n      s1 


Here  we  have  used  D  to  denote  — . 

dz 

The  boundary  conditions  are 


TJI1    =  °*  CZ1    =   °  at   z   =  -1  (3.2-17) 


T        =  0  at   z   =  s  (3.2-18) 
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At  the  interface,  z=0 


AAA 


T*1  -  Ts1  +  <M  (GZ  "  V  "  °  (3-2"19) 


A  A 


T£1  +  SeC41  +  Cl(G4  +  SeGZc-  a2A  )  -  0  (3.2-20) 


DTA1  -  BDTg1  =  oLc,  (3-2-21) 


kG*1  -  Gs1  +  h    (kGic  ~  ^    -   °  (3-2"22) 


DC 


£1  "  ^31  "  VCA1  +  vCs1  =  OV(C£c-  Csc)501         (3'2"23) 


Solving  the  system  we  obtain  a  general  equation  for  morphological 
instability 


aC.  (l-k)tanha  stanha. 
Se[G 


c       k(na  +  v/2  tanha  s)tanha„+  (a„-v/2  tanha.)tanha  s 
s         s      x-    x,         x.     s 

2.      oLtanhatanhas ,_  _  „,.,. 

-  -  GT  +  a  a  -  a(k  fcannas  +  k  tanha)       U'*-**i 

x,         s 

where         a  =  [a2  +  -  +  Y__  };2  (3.2-25) 

n   4n2 

a,  =  [a2  +  a  +  v2A]  72  (3-2-26) 

k0G.tanhas  +  k  G  tanha 

T~  k.tanhas  +  k  tanha 
9.  s 
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G<?r>(afT  v/2  tanha.)tanha  s  +  G     (na  -  v/2  tanha  s)tanha„ 

Q        _   *C   X, X, S SC    3 S   I 

c    (a  -  v/2tanha) tanha  s  +  k(na  +  v/2tanha  a) tanha. 

*  X»        S  S  3         jL 


(3.2-28) 


If  we  can  show  exchange  of  stabilities  for  this  problem,  then  for 
neutral  stability. 

GT  +  a2  A 

Se  =  (3.2-29) 

Gc 

Here  we  have  already  used  the  fact  that  the  Se  obtained  from  the  neutral 
stability  curve  will  be  the  same  as  SeQ  defined  later  in  eqn.  (4.2-2). 
In  eqn.  (3.2-29)  if  we  let  v  become  very  large  the  critical  wave  number 
amin  ^for  wnich  Se0  is  a  minimum)  also  becomes  very  large  and  we  can 
approximate  all  the  tanh  terms  to  unity.  Further  if  we  also  neglect 
solid  diffusion,  n  and  G   are  zero  and  the  equation  reduces  to  the  well 

o  (J 

known  results  obtained  by  Mullins  and  Sekerka  (1964). 

r   G8.c(a8.  "  2"  V)        aCW1"k)            ,                  2                0L 
SeQ[-^-^ S—  ♦  -l£ ]   =  t    +  a2A 2t (3.2-30) 

al   +  V(k  ~   2)  al   +  V(k   "   2J  a(kil+  k   } 

For  the  case  of  neutral  stability  this  becomes 

(1    +  a2A    )(a  /v  +  k  -  h 
Se     =  =-, ^—  (3.2-3D 

G*c(VV  "   2    ) 

It  must  be  pointed  out  that  setting  all  the  tanh  terms  to  unity  is 
equivalent  to  using  the  diffusive  length  as  the  characteristic  length. 
This  is  a  boundary  layer  approximation  not  unlike  that  used  in  the  study 
of  pipe  flow  at  high  Reynold's  numbers. 
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In  our  derivations  we  did  not  restrict  v  to  be  positive  and  hence 
eqn.  (3.2-24)  is  also  valid  for  negative  velocities,  that  is  the  fusion 
problem.  If  we  replace  v  with  -v  and  here  too  assume  that  v  is  very 
large  and  set  the  tanh  terms  to  unity  we  obtain  the  result  of  Woodruff 
(1968). 


G  (na  -V2v)        oC.  (1-k) 

0  r sc   s    c                            Jtc           l  .,    2.  oL 

SeQl ■ + J  =  1  +  a  A 

nkag  +  al  +  j  v(1-k)  nkag  +  a^  +  ^  v(1-k)  a^kx,+  k3^ 


(3.2-32) 

and  for  a  =  0 

(1  +  a2A)(nka  +  a  +  \  v(1-k)) 

Sen  =  ^—T* (3.2-33) 

Gac(l»s  -  2   v) 

To  compare  this  model  with  the  approximations  of  Mullins  and 
Sekerka  and  that  of  Woodruff,  Seomin  was  calculated  for  various  growth 
velocities  for  the  Pb-Sn  system  and  the  results  are  shown  in  Fig.  3-2. 
As  can  be  seen  the  approximations  hold  up  very  well  for  most  growth 
velocities  but  begin  to  fail  for  small  velocities.  (The  thermophysical 
data  for  the  Pb-Sn  system  were  those  of  Coriell  et  al .  (1980).)  The 
ratio  n  for  the  Pb-Sn  system  is  of  the  order  10  J  but  for  systems  with 
much  smaller  values  of  n  like  the  C-Austenite  system  (see  Clyne  and  Kurz 
(1981)  and  Wolf,  Clyne  and  Kurz  (1982))  the  approximations  begin  to  fail 
at  higher  velocities. 
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3.3.  The  Adjoint  Problem  and  Exchange  of  Stabilities 

Exchange  of  stabilities  refers  to  the  nonexistence  of  time  periodic 
infinitesimal  perturbations.  Time-dependent  infinitesimal  perturbations 
about  the  planar  state  will  generally  have  periodic  and  nonperiodic 
components,  with  o  =  a  +  ia. .  For  some  problems  it  is  possible  to  show 
that  o  is  zero  and  this  is  called  exchange  of  stabilities  (see  Iooss 
and  Joseph  (1980)  for  details). 

We  still  have  to  prove  exchange  of  stabilities  for  this  problem  but 
before  we  can  do  that  it  is  necessary  to  obtain  the  adjoint  problem.  To 
accomplish  this  we  will  define  a  column  vector  *  and  a  matrix  operator 


(3.3-1) 


L  = 


CD' 


a2) 


0 
0 


0 

Se  ,^2^2 

—  (D  -  vD-a  -  a) 

0 
0 


0 

0 

2   2 


0 
0 

0 


0     ^|(nD2-na2-a-vD) 


(3-3-2) 


and  an  inner  product     <#,*>=/   (f     T     +  C  C   )dz   +  /    (f  T     +  C  C   )dz 

i       &       Jli  JCX*  ss  ss 

-1  o 

(3.3-3) 
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where  T  denotes  the  complex  conjugate  of  T,  T*  the  adjoint  function  of 
T  and 


X  ■  (kGl<f  Gsc)/(V  V  (3-3~4) 


Then  the  domain  equations  can  be  written  as  L*.  =  0      (3.3~5) 


In  this  inner  product  the  adjoint  problem  becomes  in  the  domain 


(D^a2)^  =  0  (3.3-6) 


(D2-a2)T\  =  0  (3.3-7) 

si 


(D2  +  vD  -  a2)C*1  -  o  CJl1  (3.3-8) 


(D2+vD_a2^  =1~*  (3>3_9) 

j]  s1   n  si 


Subject  to  the  boundary  conditions 


Ti1  =  CJL1  =  °         at  Z=  ~1  (3.3-10) 


T*  =0  at  z=s  (3.3-1D 

a1 


At  the  interface  z=0 


ST    =  T  (3-3-12) 

p  1}  s1 
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nCIl  "  CIi  (3.3-13) 


kD^1  -  D^1  +  5*  (kGl0-  Gsc)  =  0         (3.3-14) 


DTIl  "D^1  +  *1  <°l"  V  =°  (3-3~15) 


DTs1  +  SeD^1  +  h    (G*  +  SeG*c"  a'A  >  =  (G^V)  ^ 

x.    S 


aC   (1-k) 

(kG,  -  G  )  SeCU  ^3.3-16) 

lO  SC 


So  far  we  have  been  unable  to  show  exchange  of  stabilities  for  this 
problem  directly.  But  Wollkind  and  Segel  (1970)  have  proved  it  when  the 
boundary  layer  approximation  is  valid  and  we  will  prove  exchange  of 
stabilities  by  performing  an  asymptotic  analysis  around  the  boundary 
layer  solution,  which  corresponds  to  n=0  and  Pe=0,  where  Pe  is  the 
Peclet  number  given  by  D./Jlv. 

JO 


T   =  T°°  ♦  nT1°  ♦  PeT01   ♦  (3  3-17) 


a°°  +  no10  +  Pea01  +  (3-3-18) 


The  other  variables  are  expanded  similarly. 

If  we  use  L   to  designate  the  operator  L  in  (3.3~2)  when  n  and  Pe  are 

zero,  then  the  linear  perturbed  systems  become 


l00  ;;°  =  f10     &       l°°  ;°1 .  f°1      0.3-19)  &  u.3-20) 
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where 


.10 


0 
0 

10:00 


a  C 


8,1 


.01 


0 
0 

01200 

a  C  , 
s1 


(3.3-21)  &  (3.3-22) 


With  boundary  conditions 


:io     :io 


at  z  =  -  • 


(3.3-23) 


?°  =  0 
s1 


at  z=°> 


(3.3-24) 


at  z=0  the  equations  are 


D00,^10  "10,   JO 

B  (^  ,  c1  )  =  h      & 


B00(;;°,  ^>  =h01   (3.3-25)  &  (3.3-26) 


where  B   is  the  boundary  operator  defined  by  eqns.  (3.2-19)  -  (3.2.23) 
when  n  and  Pe  are  zero  and 


10 


poo     10_ 

C1    (GZ 

Koo, 


c> 


10*00 
a     LC1 


Z°™i°  ■■  O 


■i  v    ac 

10^00. 
'Hc( 


,     "00  00n10,,    ,     "00 

'     ,(1-k)51      +   o      Cic(1-k)^1 


(3.3-27) 
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h  will  have  a  similiar  form  with  01  superscripts  replacing  the  10.  It 
we  let  L  and  ♦  represent  the  adjoint  operator  and  adjoint 
function  of  L  and  ♦  respectively,  we  can  use  the  solvability 
condition  on  the  above  system. 


<;00*f  L00;i0>  m   ^00^  f10>  (3>3_28) 


oo*~oo*  "in 
<LUU  #"°  ,  *11U>  =  0  (3.3-29) 


Subtracting  we  get 


T,:oo*  "oo*  jo, 

J(*1   ,  C1   »  h  ) 


=  <*°0*,  f01>       (3.3-30) 


z=0     1 


where  the  lefthand  side  of  eqn  (3.3~30)  is  the  bilinear  concomitant 
evaluated  at  z=0. 


Similiarly  J(#   ,5   ,  h  ) | z=0  ■  <*1   »  f  >         (3.3-31 ) 


~00*     "00 
We  note  that  *    and  c,   are  real  as  they  correspond  to  a  state  where 

exchange  of  stabilities  has  been  proved,  while  h  ,  h  ,  f   and  f   are 

also  made  up  of  real  quantities  except  for  o   and  o  .  Hence  we 

conclude   from   eqn   (3.3-30)   that  a   is   real   and  from   (3.3-31 ) 

that  a   is  real;  i.e.  exchange  of  stabilities  holds  for  the  generalized 

morphological   instability   problem  upto  0(n)  and  0(Pe).    So  we  are 

justified  in  using  the  neutral  stability  curve  (3-2-29)  to  calculate  SeQ 

at  least  up  to  such  order  even  though  we  often  extend  it  further.  Even 

when  the  boundary  layer  approximation  is  valid,  exchange  of  stabilities 
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does  not  hold  for  all  growth  conditions  (see  Coriell  and  Sekerka  (1983)) 
and  care  must  be  exercised  when  such  extensions  are  made. 

3.4.  Finite  Containers  and  the  Most  Dangerous  Wavenumber 

Appealing  to  the  proof  of  Section  3-3,  henceforth  we  shall  only 
consider  steady  state  solutions.  A  typical  Se  versus  a  diagram  is 
shown  in  Fig.  3-3  and  Seomln  is  the  minimum  value  of  SeQ  and  the 
wavenumber  at  which  this  occurs  is  amin.  If  we  can  maintain  the  growth 
conditions  such  that  Se  <  SeQmin  we  can  at  least  say  that  the  planar 
interface  is  stable  to  infinitesimal  perturbations. 

As  Seomin  is  the  least  value  of  Se  at  which  the  planar  solution 
loses  the  stability,  amin  is  the  wavenumber  of  the  disturbance  which  is 
most  likely  to  occur.  Hence  this  wavenumber  is  commonly  regarded  as  the 
most  dangerous  and  is  the  wavenumber  at  which  morphological  instability 
is  usually  studied.  In  the  derivations  and  discussions  that  follow  this 
is  the  wavenumber  employed  and  it  becomes  our  "operating  wavenumber." 

The  wavenumber  is  a  2ir  multiple  of  the  reciprocal  of  the 
wavelength.  If  the  domain  being  considered  is  regarded  as  being  finite 
with  periodic  lateral  boundary  conditions,  the  wavelength  R  (or  L 
depending  on  the  wave  pattern  used)  in  its  dimensionless  form  is  now  the 
aspect  ratio,  the  ratio  of  the  ampoule  radius  to  the  melt  depth.  For 
this  situation  SeQ  takes  on  different  values  depending  on  the  number  of 
cells  formed  on  the  crystal  surface  as  shown  in  Fig.  3-4  (see  also 
Rosenblat,  Homsy  and  Davis  (1982)).  For  each  value  of  R  there  is  a 
fixed  number  of  cells  which  is  the  pattern  that  is  most  easily 
disturbed,  except  at  certain  values  of  R  where  two  different  patterns 
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amln 
WAVENUMBER  a 


Figure  3-3.   Se  vs.  wavenumber  diagram 
o  ° 


n=51 


n  =  52 


Figure  3-4.   Se  curves  for  various  aspect  ratios 


37 


are  equally  dangerous.  These  "horizontal  multiple  points"  are  very 
important  and  will  be  discussed  later.  At  other  values  of  R  the  most 
dangerous  number  of  cells  is  denoted  by  N  and  the  corresponding 
wavenumber  of  each  cell  and  the  value  of  SeQ  are  denoted  as  aN  and  SeQN 
respectively.  For  the  analysis  in  Chapter  6  where  the  effect  of  a 
finite  container  on  morphological  instability  is  considered,  R  is  chosen 
such  that  SeQN  is  as  close  as  possible  to  Seom^n  so  that  the  worst 
possible  case  can  be  examined. 


CHAPTER  H 
SUBCRITICAL  BIFURCATIONS 


4.1 .  Theory 


The  linear  stability  analysis  will  only  give  the  onset  conditions 
for  morphological  instability.  Nonlinear  calculations  are  necessary  to 
determine  the  behavior  beyond  this  point.  Ideally  numerical 
calculations  should  give  the  most  amount  of  information  by  being 
applicable  for  small  and  large  deformations,  but  as  can  be  seen  from  the 
work  of  Ungar  and  Brown  these  calculations  are  very  expensive  to  carry 
out  and  they  were  forced  to  simplify  the  nonlinear  problem  and  perform 
calculations  for  very  few  experimental  conditions.  McFadden  et  al . 
(1987),  even  though  they  did  not  attempt  calculations  for  larger 
deformations,  were  faced  with  the  same  restrictions. 

This  then  is  the  case  for  weakly-nonlinear  methods.  They  are 
generally  valid  only  in  a  small  region  very  close  to  the  onset 
conditions  but  they  can  be  used  to  predict  the  shape  of  the  nonlinear 
curve  for  larger  deformations  and  for  several  applications  this 
information  is  sufficient.  More  importantly  due  to  the  analytical 
nature  of  the  techniques  they  can  be  used  to  predict  the  weakly 
nonlinear  behavior  for  all  experimental  conditions.  A  case  in  point  is 
the  work  of  McFadden  et  al . ,  most  of  whose  predictions  could  have  been 
obtained  more  cheaply  for  all  parameter  values  from  weakly-nonlinear 
theories. 
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Probably  the  most  useful  information  generated  by  these  theories  is 
the  subcritical  behavior  of  the  nonlinear  curve.  Some  typical 
bifurcation  diagrams  are  shown  in  Figs.  4-1  to  4-4.  The  e  =  0  axis 
corresponds  to  the  planar  solution  and  initially  for  small  values  of  Se 
the  planar  solution  is  stable  and  usually  the  only  possible  solution. 
For  Se  £  SeQN  the  planar  solution  becomes  unstable  and  a  nonplanar 
solution  "bifurcates"  from  the  planar  one.  Figure  4-1  shows  a  symmetric 
bifurcation  diagram  where  nonplanar  solutions  do  not  exist  for  Se  < 
SeoN*  For  Se  >  SeoN'  even  an  infinitesimal  perturbation  will  make  the 
solution  jump  from  the  unstable  planar  solution  to  the  stable  planar  one 
while  for  Se  <  SeoN  the  planar  solution  is  stable  to  all 
perturbations.     This   behavior   is   referred   to   as   supercritical 

bifurcation  and  for  these  curves  Se1  =  0,  Se2  >  0  (where  Se1  =  dSe/de 

2     2 
and  Se2  =  d  Se/de  at  e  =  0).   Figure  4-2  is  nonsymmetric  and  as  can  be 

seen  stable  and  unstable  nonplanar  solutions  exist  for  Se  <  Se _«,  making 

the  planar  solution  stable  to  infinitesimal  perturbations  in  this 

region,  while  a  large  perturbation  can  make  it  jump  to  the  stable 

nonplanar  branch.   This  is  called  a  subcritical  bifurcation  diagram  and 

is  characterized  by  Se1  *  0.   In  this  situation  it  is  obvious  that 

growing  the  crystal  at  Se  <  SeQN  is  no  guarantee  of  avoiding 

morphological  instability. 

Figures  4-1  and  4-2  display  the  behavior  usually  seen  in  most 

problems  of  hydrodynamic  instability.    Morphological  instability  is 

unusual  in  having  nonlinear  curves  shown  by  Figs.  4-3  and  4-4  as  well. 

These  curves  have  been  labelled  "backward  bending"  to  distinguish  them 

from  the  usual  "forward  bending"  curves.   (Actually  they  are  Janus-like 

in  appearance  bending  backwards  and  forwards.)   From  the  point  of  view 
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Figure  4-1.  Forward  bending,  locally 
symmetric,  supercritical 
bifurcation  diagram 


Figure  4-2. 


Forward  bending, 
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critical  bifur- 
cation diagram 
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Figure  4-3.   Backward  bending,  locally 
symmetric,  subcritical 
bifurcation  diagram 


Figure  4-4. 


Backward  bending, 
unsymmetric,  sub- 
critical  bifurca- 
tion diagram 
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of  the  crystal  grower  this  is  unfortunate  as  their  subcritical  nature 
increases  the  occurrence  of  subcritical  bifurcation.  Figure  4-3  is  the 
symmetric  case  with  Se-j  =  0  and  Se2  <  0  and  Fig.  4-4  the  nonsymmetric 
one  with  Se  *  0. 

The  symmetry  or  nonsymmetry  of  the  bifurcation  diagram  in 
hydrodynamic  stability  is  dependent  on  the  cell  pattern  (see  Joseph 
(1976)).  For  morphological  instability  it  will  be  shown  that  two 
dimensional  rolls,  rectangular  cells  and  cross  rolls  produce  locally- 
symmetric  bifurcation  diagrams  while  cylindrical  rolls  and  hexagonal 
cells  produce  nonsymmetric  diagrams.  It  will  also  be  shown  that  both 
backward  bending  and  forward  bending  will  occur  for  all  cellular 
patterns  depending  on  the  experimental  parameters  used.  Hence 
bifurcation  can  be  subcritical  or  supercritical  for  two  dimensional 
rolls,  rectangular  cells  and  cross  rolls  but  for  hexagonal  cells  and 
cylindrical  cells  bifurcation  is  always  subcritical. 

4.2.  The  Second  Order  Problem 

Here  we  begin  our  weakly  nonlinear  analysis  and  in  this  section  we 
will  calculate  Se1 ,  the  first  derivative  of  Se  with  respect  to  e. 
Considering  the  neutrally  stable  nonplanar  solution  near  the  bifurcation 
point  Se  ,  we  will  expand  the  variables  around  the  planar  steady  state 
solution. 


T  =  T   +  eT  +  e  T   +  £  T  +  (4.2-1) 

I         He     Jto     11      12 


Se  =  Se  +  eSe,  +  e2Se_  +.  (4.2-2) 

o      1       2 
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where  e  -  <♦  -  *  ,  +  -  *  >  f2  (4.2-3) 

We  have  already  obtained  the  linear  perturbed  solution  in  Section 
3.2.  Substituting  these  expansions  into  the  steady  state  versions  of 
eqns.  (3.1-15)  -  (3.1-24)  and  collecting  the  terms  of  order  e  we  get 
the  second  order  perturbed  problem. 

If  the  first  order  perturbed  variables  are  written,  following  eqn. 
(3.2-12),  as 


T&o(x,y,z,Seo)  =  T^  (z.Se^  (x,y)  (4.2-4) 


Then  the  solution  to  the  second  order  problem  becomes 


T*1  "^  TA1n(z'SV  W»n(x'y)  (4'2-5) 


Substituting  (4.2-4)  and  (4.2-5)  into  the  second  order  problem  and 
taking  Fourier  transforms  horizontally,  we  get 


l*u    =0  (4.2-6) 


and  TA1 1  =  C£1 1  =  °        at  z  =  ~1  (4.2-7) 


Tal1  =0  at  z  =s  (4.2-8) 


B(*n,  cn)  =  hn       at  z=0  (4.2-9) 
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where  L  and  B  are  the  same  as  that  used  in   (3.3-2)  and   (3«3-25)   but  with 

o=0  and  Se  =  Se     and 
o 


11 


-C      (DT  -   DT        )I 

6or    ioi         sor  11 


CmGflJ 


^01(DT*01    +  SeoDW    -   2  ^^VOI^II    -   Se1(C£01 

"a%i(T£or  ^aoi^n  +  WW  aW1^ 
[-WkDC*oi  -DW  -  I  Y*m{Mt*  -  am/  n)]ln 


"a  Wcioi  -  ^.oi^n  + 


501(CJt01  "  nCa01)Il2 


(4.2-10) 


■11 


L2   L1 


/  /   4^(x,y)dxdy 


o  o 


L2L1 


J  /  <()1(x,y)dxdy 


o  o 


(4.2-11) 


L2  S    9<J»       3« 


o   o 


3x 


3y 


'12 


L2L1 


/  J  4  (x.y)dxdy 


o  o 


(4.2-12) 


It  can  easily  be  seen  that  for  two  dimensional  rolls,  cross  rolls 

and  rectangular  cells  that  In  =  I12  =  0.   Hence  from  the  solvability 

condition  for  these  patterns  Se-,  will  be  zero;  that  is,  the  bifurcation 

will  be  locally  symmetric.   But  for  hexagonal  cells  and  cylindrical 
rolls  1^  and  I12  will  be  non-zero  and  hence  Se1  will  also  be  nonzero 

and  we  have  nonsymmetric  bifurcation  and  the  existence  of  subcritical 
instabilities . 
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We  will  analyze  the  case  of  Se  *  0  by  considering  hexagonal  cells 
as  an  example,  but  the  results  obtained  are  applicable  to  cylindrical 
rolls  as  well.  It  should  be  noted  here  that  all  the  horizontal  cell 
patterns  we  have  been  using  are  possible  solutions  to  Helmholz's 
equation 


4+4+*2*=0  (4'2"13) 

dx         9y 

and  the  solutions  always  come  in  pairs.    For  hexagonal  cells  the 
complementary  solution  to  <\>     given  in  (3.2-10)  is  ty 


♦1  -  2  Cos  7|j-x  Sin  -|£y  -  Sin  |jy  (4.2-14) 


So  the  general  form  of  eqn.  (5.4)  will  be 


T*o  =  W*i  +>V  (4-2"15) 


To  determine  p  the  procedure  outlined  by  Joseph  (1976)  will  be 
used.  We  will  proceed  the  same  way  as  above  but  using  (4.2-15)  instead 
of  (4.2-4)  and  multiplying  by  <p.  and  integrating  horizontally.  This 
will  give  the  same  set  of  eqns.    (4.2-6)   -   (4.2-10)   but  with 

)      o 

3~lT73L 

J      J      (♦1   +  pipj^dxdy 


3L  /3L 

J     J      (^    +  P^)   ^dxdy 

o     o  , ,      .2 


rl,  -1T73L-  ■-(,-")  <"-2-16) 


O       0 
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3L  /3L         3<j>  3tK  3*  3i|i 

T  o     o a     ,.      2. 

I12  -  3L-73L =  —  Cl-p   ) 

J      J      (♦1    +  P^^dxdy 
o     o 

(4.2-17) 


which  can  be  used  to  calculate  Se,  .  We  then  repeat  the  process 
with  \\>.  and  equate  the  two  Se  's  obtained.  This  will  result  in  a  cubic 
equation  for  p. 


p3  -  3p  -  0  (4.2-18) 


and  Se   =  Se. (1-p2)  (4.2-19) 


Se,  G  (kG„   -  G  )   (a.  +  v/2tanha.)B   (na  -  v/2tanha  s) 

1   C  =  -      JO     SC   jr  I l_    +     3 S -I  j 

C01    ■   o  (k  +  g)2       tanha^  ntanha^       11 


(1+nB)tanha  s 

(na  ;  v/2tanha  s)  W  +  h1  ^^^ 

s         s 


»<Gt  ~   V  r  aI11  +  aiGl   +  SeoG£c  ~  a2A)(1/S-D 

h1  =  (B+A)     '■tanha  +       (G.-G  )tanha       "  11 

l     s 


(1 +A)tanhas 

(B+A)       "12  J   "2 


I„  }  ♦  h0  (4.2-21) 


Se  I,,   a(B-1)(kG.   -  G  ) 

h  =   o11[- - i2 sc__(     +G  /}         (4.2-22) 

2   (k+B)  L  (B+A)tanha  4c    sc 
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where  B  =  (a„  -  v/2tanha„ )tanha  s/(na  +  v/2tanha  s)tanha.    (4.2-23) 
x.  *      s     s  3       x, 


and    A  =  tanhas/tanha  (4.2-24) 

There  are  three  solutions  0,  /3  and  -/3  for  p,  but  it  can  be  seen 
that  Se-|D  for  /3  and  -/3  coincide.  It  is  also  obvious  that  while  any 
two  of  the  solutions  are  independent  the  third  is  a  linear  combination 
of  the  other  two.  Here  is  an  instance  where  the  problem  exhibits  a 
multiplicity  of  solutions  for  the  same  eigenvalue,  and  could  cause 
secondary  bifurcations  further  along  the  bifurcation  curve. 

The  next  obvious  question  is  to  ask  if  there  is  any  point  at  which 
Sen  in  (4.2-20)  goes  to  zero  and  hence  causing  Se1p  to  go  to  zero. 
Calculations  done  for  the  C-Austenite  (i.e.  steel)  system  did  give  such 
a  curve  (see  Fig.  4-5)  though  for  most  growth  conditions  this  curve  lies 
far  away  from  the  critical  (or  most  dangerous)  wavenumber  amin, 
intersecting  it  only  at  high  velocities. 

4.3.  The  Third  Order  Problem 

In  Section  4.2,  it  was  shown  that  for  2-D  rolls,  cross  rolls  and 
rectangular  cells  Se«  was  zero,  which  implied  a  symmetric  bifurcation 
diagram.  But  to  learn  more  about  the  nature  of  this  diagram  it  is 
necessary  to  go  to  the  next  order.  If  we  repeat  the  procedure  described 
for  the  second  order  problem  for  the  terms  of  order  £-*  we  will  obtain 
the  third  order  problem. 


L*21  =  0  (4.3-1) 
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Figure  4-5.      Growth  velocity  vs.   wavenumber  chart   for 
C-Austenite   system  for  hexagonal  cells 
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T  =    C  =    0 


at  z=  -1 


(4.3-2) 


Ts21    "  ° 


at  z=s 


(4.3-3) 


B(#21    ,   t     )   =h21 


at  z=0 


(4.3-4) 


21 


a2  2 

-  -  C      (T  -   T        )I 

2  ^or  401       sor  21 

a2     2 


v   2 


1    2,3 


■V   I" ^01CTZ01    +   SeoW    +   2^01Seo   DC*01    +   6V   'o>oGJ 


"    (GZ0t   ^01G&c)Se2  "     A4(3/2a2l22"   I23) 


a2     2 


'^iHt  coi(kC*oi  -  W  +  f^i(k  DC*oi  "  DCsoi/rl) 

+  ^01(kGJlc-Gsc/ri2)J 


2 
I  ?01V(DCZ01 


nDCs01    U21    +  <VDCil01    "  nDC301)I22 


L2     L1 


J     /     *J    dxdy 


o     o 


"21 


L2L1 


|     /  4>2dxdy 


o     o 


(4.3-5) 


(4.3-6) 


L2  L1        2     3*1      2        3<{. 


/  J    ♦fK-k   +  (4v}  ]dxdy 


o     o 


"22 


L2      L1 


/       /  ♦fdxdy 


o       o 


(4.3-7) 
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L\    L2     3<J>         3*  3*     3$     22<|>         32<(,        3<(, 

:23  =  TTTT 

o     o 


(4.3-8) 


Once  again  using  the  solvability  condition  we  obtain 


Se  2  v(kG.   -G     )(1+B)tanha  s 

J-  9  ■  „    «  „    .  (1      -  |-  !     )s«  «o    -  »        (4.3-9) 

Cn,  (k+B)    (na  +1/2vtanha  s) 

01  S3 

where 


f32               w                       rvGic       (kGjic"G3c)(a2.+  v/2tanha   )/tanha 
h1    "  ^1*   I22+I23^-   la*Q"   t-T  "  6(k+B) * ^ 


(4.3-10) 


.  *21    V(G&C-GSC)       rk(V   V/2   ^V 

2"         2  (k+B)2       I  tanha^ 


(na  -v/2  tanha  s)B  Se  v(kG.   -G     /n) 

s  a  l    „  T  O  SoC      SC  /  h    o    1  1  i 

2— J  Seo  +  1^^  6TkTB) (iK3-11) 

n  tanha  s 

s 


and  B  =•   (a .-1/2  vtanha. ) tanha  s/(na  +   1/2  vtanha  s)tanhan  (4.3-12) 

x,  x,  s  s  s  8, 
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4.4.  Calculations  and  Comparisons 

Using  the  results  of  the  previous  section,  Se2  was  calculated  for 
various  values  of  the  experimental  parameters  and  the  results  are  shown 
in  Figs.  4-6  to  4-9.  Figure  4-6  was  drawn  in  order  to  check  the 
derivation  with  that  of  Ungar  and  Brown  (1984a)  and  shows  calculations 
done  for  2-D  rolls  with  Dg  =  0  and  8  =  1  but  is  a  more  complete 
calculation  for  the  experimental  parameter  ranges  involved  than 
theirs.  The  calculations  agree  with  those  of  Ungar  and  Brown  but  it  is 
an  unexpected  result  nevertheless,  showing  multiple  regions  of 
subcritical  and  supercritical  instability.  It  also  seemed  to  contradict 
the  earlier  calculations  of  Wollkind  and  Segel  (1970)  who  did  not  see 
any  supercritical  instability  but  this  was  resolved  in  the  paper  of 
Wollkind  and  Wang  (1988)  and  hence  Fig.  4-6  agrees  with  their 
calculations  as  well.  As  argued  in  Section  3.4  if  we  treat  the  cl, 
curve  as  the  experimental  "operating  line,"  the  bifurcation  is  mostly 
supercritical  which  also  seems  to  contradict  the  experimental  results  of 
de  Cheveigne  et  al.  (1986)  who  observed  only  subcritical  bifurcation. 

In  Fig.  4-7  calculations  were  done  for  more  realistic  values  of  DG 
and  0  and  the  results  are  significantly  different  with  only  one  region 
of  supercritically  and  another  of  subcriticality.  Ignoring  solid 
diffusion  introduces  a  singularity  to  the  problem  and  this  accounts  for 
the  distortions  observed  in  Fig.  4-6.  In  Fig.  4-7  if  we  move  along  the 
operating  line  for  a  fixed  liquid  temperature  gradient,  initially  for 
small  growth  velocities  the  bifurcation  is  supercritical,  until  at  a 
critical  velocity  a  transition  point  is  reached  and  the  bifurcation 
becomes  subcritical.    Hence  for  every  imposed  liquid  temperature 
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Figure  4-6.   Velocity  vs.  wavenumber  chart  for  the  one-sided 
approximations  of  Ungar  and  Brown  for  2-D  rolls 
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gradient  there  will  be  a  critical  growth  velocity  below  which  the 

bifurcation  for  2-D  roll  disturbances  is  always  supercritical.   The 

really  surprising  result  here  is  that  when  these  transition  points  for 

various  values  of  the  gradient  are  joined,  we  get  a  straight  line 

through  the  origin.   We  now  have  a  clearly  demarcated  supercritical 

upper-triangular  and  subcritical  lower-triangular  one.  It  is  well  known 

that  the  onset  condition  SeQ  does  not  change  much  in  the  neighborhood  of 

^in  (see  for  sample  Coriell,  McFadden  and  Sekerka  (1985))  and  so  a_. 

min 

is  more  of  an  interval  than  a  unique  point.  It  can  also  be  seen  from 
Fig.  4-7  that  the  amin  curve  practically  hugs  the  line  of  transition 
from  sub  to  supercriticality  and  if  we  impose  an  interval  for  a  .  it 
would  straddle  the  transition  line.  Thus  for  2-D  rolls  it  is  unlikely 
that  we  will  ever  see  a  sharp  transition  from  subcritical  to 
supercritical  bifurcation  in  experiments.  More  likely  we  will  observe 
subcritical  behavior  throughout  as  reported  by  de  Cheveigne  et  al . 

Proceeding  to  the  three  dimensional  patterns,  we  obtained  almost 
identical  results  for  square  cells  and  cross  rolls.  Figure  4-8  is  for 
square  cells  and  we  see  quite  a  change  with  the  supercritical  region 
acquiring  a  characteristic  balloon  shape  and  having  a  sharp  transition 
to  subcritical  bifurcation  along  the  amin  curve.  But  here  too  if  these 
points  of  transition  are  joined,  a  straight  line  is  the  result, 
demarcating  a  supercritical  upper-triangular  operating  region  and  a 
subcritical  lower  triangular  one.  Unlike  the  case  of  2-D  rolls  these 
should  be  visible  to  the  experimentalist.  So  in  order  to  avoid  cross 
rolls  and  square  cellular  instabilities  not  only  should  the  crystal  be 
grown  when  Se  <  Seomin,  but  one  should  do  so  in  the  upper  triangular 
region.  Figure  4-9  demonstrates  the  universality  of  our  result  in  being 
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applicable  for  the  melting  problem  as  well  and  repeating  the  derivation 
for  this  case  separately  as  was  done  by  Wollkind  and  Raissi  (1974)  is 
unnecessary. 

Finally,  even  though  for  hexagons  Se1  was  usually  non-zero  (and 
hence  Se2  cannot  easily  be  calculated),  in  Section  4.2  we  saw  that  there 
were  points  at  which  Se1  did  go  to  zero.  If  we  attempt  to  evaluate  Se2 
for  hexagons  at  these  points  I21 ,  I22  and  I2«  become 


J21  -  nr  (1  +  P2)  (4.4-1) 


I22  =fa2(1+p2)  (4.4-2) 


I23  =  -  |  a  (1+p2)  (4.4-3) 


and  as  can  be  expected  Se2  will  be  in  the  form 


Se2p  =  Se2(1  +  p2)  (4.4-4) 


where  Se2  is  that  corresponding  to  p  =  0.  The  other  possible  value  is 
when  p  is  /3  or  -/3.  Here  Se2p  can  be  calculated  from  (4.3-9)  -  (4.3- 
12)  and  (4.4-1)  -  (4.4-3).  As  mentioned  in  Section  4.2  these  points  are 
usually  far  away  from  amin  and  in  Fig.  4-5  we  found  that  along  this 
curve  Se2  is  positive  at  low  values  of  "a."  As  "a"  increases,  at  a 
point  above  a=amin,  Se2  becomes  negative.  Hence  along  this  line  the 
bifurcation  diagram  is  forward  bending  at  low  "a's"  (including  a^.  )  and 
becomes  backward  bending  at  high  values  of  "a."   Even  though  this  is 
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true  only  along  the  Se.,=0  line,  by  analogy  with  other  cell  patterns  we 
conjecture  that  this  is  valid  when  Se  *0  as  well.  In  other  words  we 
expect  the  bifurcation  diagram  to  be  forward  bending  along  the  a_.  line 
and  below  (as  shown  in  Figs.  4-10  and  4-11)  but  a  transition  to  backward 
bending  along  the  amin  line  could  occur  at  high  velocities.  Far  above 
the  amin  line  the  bifurcation  diagram  should  be  backward  bending  (see 
Figs.  4-12  and  4-13). 

To  confirm  these  conjectures  we  turned  to  the  nonlinear 
calculations  for  hexagons  of  McFadden  et  al.  (1987)  but  unfortunately 
they  were  unable  to  complete  the  bifurcation  diagram  as  their  attempts 
to  compute  the  curve  for  e<0  failed.  Also  their  calculations  were  done 
only  for  the  case  of  p=0  and  not  for  p=±J?>.  But  they  did  confirm  the 
existence  of  subcritical  instability. 

To  sum  up,  it  has  been  shown  that  the  Mullins  and  Sekerka  and  the 
Woodruff  models  of  morphological  instability  are  of  limited  validity. 
The  uniform  formulation  is  the  more  exact  representation  of  the  problem 
and 

*  it  is  applicable  for  all  growth  velocities  and  not  just  the 
relatively  rapid  solidification  and  fusion  regions  and  provides  a 
single  formulation  from  which  all  the  different  models  for 
various  growth  conditions  can  be  obtained  as  limiting  cases, 
eliminating  the  duplication  of  derivations  for  different  cases; 

*  it  incorporates  the  whole  crystal  and  melt  regions  into  the 
problem  and  not  just  a  boundary  layer  region  adjoining  the 
interface  facilitating  the  study  of  various  domain  effects  like 
convection  on  morphological  instability; 
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Figure  4-13.   Bifurcation  diagram 
for  hexagonal  cells 
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*  it  avoids  the  incorrect  predictions  of  subcritical  bifurcation 
regions  because  of  the  singularities  inherent  in  previous  models. 

*  The  principle  of  exchange  of  stabilities  has  been  shown  to  be 
applicable  to  this  model  as  well  even  though  only  in  an 
asymptotic  sense. 

When  the  weakly  nonlinear  technique  of  Malkus  and  Veronis  (1959)  is 
applied  to  this  problem  in  a  systematic  way,  it  resulted  in  important 
information  about  the  shape  of  the  bifurcation  diagram  for  various 
growth  conditions.  Some  of  these  results  are  similar  to  those  obtained 
for  Marangoni  instability  (see  Joseph  (1976))  which  leads  us  to  assert 
that  these  results  are  valid  for  all  hydrodynamic  instability  problems 
in  which  the  nonlinearity  lies  only  on  the  boundary. 

*  If  rectangular  cells,  cross  rolls  or  two  dimensional  rolls  are 
the  horizontal  cell  patterns  then  the  bifurcation  diagram  will 
always  be  locally  symmetric.  For  hexagonal  cells  or  cylindrical 
rolls  they  are  generally  non-symmetric  and  hence  the  bifurcation 
is  subcritical. 

*  The  problem  can  display  a  multiplicity  of  solutions  for  the  same 
eigenvalue.  Specifically  for  hexagons  there  are  two  possible 
solutions. 

Considering  the  morphological  instability  problem  in  particular  the 
following  were  shown: 

*  for  two-dimensional  rolls  there  are  two  operating  regions,  one 
subcritical  (backward  bending)  and  the  other  supercritical 
(forward  bending),  but  since  the  demarcation  is  not  sharp  its 
probable  that  only  subscritical  bifurcation  will  be  observed 
experimentally. 
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*  For  rectangular  cells  the  forward  bending  region  has  a 
characteristic  balloon-like  shape  and  here  too  there  is  a 
straight  line  dividing  the  operating  region  into  subcritical  and 
supercritical  zones,  but  here  the  transition  is  sharp  and  hence 
probably  observable  experimentally. 

*  For  hexagonal  cells  and  cylindrical  rolls  the  bifurcation  diagram 
shows  both  backward  and  forward  bending  behavior  but  the  exact 
regions  of  each  can  only  be  conjectured. 


CHAPTER  5 
COMPARISONS  WITH  RAYLEIGH-MARANGONI  CONVECTION 


5.1.  Rayleigh-Marangoni  Convection  in  Brief 

When  a  horizontal  layer  of  quiescent  fluid  is  heated  from  below,  on 
account  of  thermal  expansion,  the  fluid  at  the  bottom  will  be  lighter 
than  the  fluid  at  the  top.  This  unstable  arrangement  is  maintained  by 
the  viscosity  of  the  fluid  which  inhibits  any  flow  and  suppresses 
disturbances  such  that  there  will  be  a  stable  conduction  profile  in  the 
fluid.  But  when  the  adverse  temperature  gradient  exceeds  a  critical 
value,  the  viscous  force  is  overcome  and  there  will  be  cellular 
convection.  This  instability  is  known  as  Rayleigh-Benard  convective 
instability. 

There  are  several  variations  of  this  problem.  Instead  of  an  ad- 
verse temperature  gradient,  there  could  be  an  adverse  solute  concentra- 
tion gradient  in  the  fluid  causing  once  again  an  unstable  top-heavy 
arrangement.  The  convective  instability  arising  from  this  is  known  as 
solutal  convection  or  the  solutal  Rayleigh  problem.  Another  way  to 
cause  convection  is  to  have  a  very  thin  fluid  layer  heated  from  below, 
but  the  top  surface  of  the  fluid  instead  of  being  kept  at  a  fixed  lower 
temperature,  is  allowed  to  remain  free.  Here  the  thinness  of  the  fluid 
layer  makes  buoyancy  effects  negligible  but  convection  will  be  caused  by 
surface  tension  variation  on  the  free  surface.  This  is  known  as  Maran- 
goni  convection.   Finally,  there  could  be  combinations  of  the  above 
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three  types  of  convective  instability.  When  convection  is  caused  by 
thermal  and  concentration  gradients  it  is  known  as  double-diffusive 
convection.   Combination  of  either  thermal  or  solutal  convection  with 
surface  tension  driven  flow  is  the  Rayleigh-Marangoni  problem. 

In  addition  to  these  there  are  several  other  combinations  possible 
like  Rayleigh-Benard  convection  with  rotation  or  with  a  magnetic  field 
but  for  purposes  of  comparison  with  morphological  instability  it  will  be 
seen  that  the  three  causes  for  convection  mentioned  above  are  the  most 
relevant. 

In  the  discussion  to  follow  it  is  desirable  to  consider  the  most 
general  form  of  this  problem.  Despite  there  being  several  interesting 
features  in  the  problem  of  double-diffusive  convection,  the  causes  for 
convection  there,  the  temperature  gradient  and  the  solute  concentration 
gradient  are  both  bery  similar  and  it  is  sufficient  to  look  at  the 
effect  of  one  gradient.  The  manner  in  which  the  surface  tension  effects 
the  problem  is  very  different  from  the  buoyancy  effects  and  a  general 
formulation  should  include  both.  Accordingly  we  will  examine  the  Ray- 
leigh-Marangoni problem  with  thermal  convection.  In  the  following 
sections  several  other  reasons  for  looking  at  this  problem  will  become 
apparent. 

The  equations  for  the  Rayleigh-Marangoni  problem  are  given  by  Sarma 
(1987)  and  we  will  reproduce  them  here  and  refer  the  interested  reader 
to  his  paper  for  details.  The  steady  state  dimensionless  Boussinesq 
equations  in  the  domain  are 

V-V  =  0  (5.1-1) 
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2         Ra 
V  V  "  VP  +  p^  Tg  -  V-VV  =  0  (5.1-2) 


V2T  -  PrV'VT  =  0  (5.1-3) 


where  p  is  the  modified  pressure,  g  the  acceleration  due  to  gravity,  Ra 
the  Rayleigh  number  and  Pr  the  Prandtl  number.  The  boundary  conditions 
at  the  bottom  of  the  fluid  layer  are 

at  z  =  0,  T  =  T  and  w  =  0,  8w/3z  =  0  (5.1-4) 

where  w  is  the  vertical  component  of  velocity. 

The  boundary  conditions  at  the  top  are  more  complicated.  Not  only 
will  there  be  surface  tension  variation  across  the  surface  but  the 
surface  is  also  free  to  deflect  like  the  liquid-solid  interface  in 
crystal  growth. 

at  z  -  1  +  c,         VT«n  =  -  BiT  (5.1-5) 

v*n  =  0  (5.1-6) 


Bopn  +  Crn«[VV  +  VVT]  =  MaV  T  +    f/n  (5.1-7) 

H 


The  dimensionless  quantities  are  Bi  the  Biot  number,  Bo  the  Bond  number, 
Cr  the  Crispation  number,  Ma  the  Marangoni  number  and  H  the  surface 
curvature  (see  Scriven  and  Sternling  (1964)  for  details). 

Initially  there  will  be  a  quiescent,  linear,  stable,  conducting 
solution  to  the  problem  with  V  =0.   At  a  critical  value  of  the  charac- 
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teristic  parameter  (Ma  or  Ra)  this  conducting  solution  becomes  unstable 
to  infinitesimal  perturbations  and  we  have  a  convective  solution. 
Performing  a  linear  stability  analysis  about  the  conduction  state, 
separating  variables  and  doing  considerable  manipulations  we  get  in  the 
domain 


2    2  3       2 
(D  -  a  )Je  =  -  a  Rae  (5.1-8) 


2    2  3       2 
(D  -  a  )"w  =  -  a  Raw  (5.1-9) 


where  "a"  is  the  wavenumber,  w  the  Fourier  coefficient  of  the  vertical 
component  of  velocity  and  e  the  Fourier  coefficient  of  the  temperature. 

At  the  boundary  at  z  =  0, 

w  =  Dw  =  6  =  0  (5.1-10) 

At  z  =  1,  w  =  0  (5.1-11) 


2     2 
BiD  w  =  a  MaD6  (5.1-12) 


BiCr(D3w  -  3a2Dw)  =  a2(Bo  +  a2)(D6  +  Bi9)        (5.1-13) 


When  the  density  variation  with  temperature  is  negligible  or  in  the 
absence  of  gravity,  then  Ra  =  0  and  we  have  the  pure  Marangoni  problem 
with  all  the  nonli near i ties  only  in  the  boundary.  Even  with  this  effect 
in  the  boundary  the  important  result  of  bifurcation,  namely  the  fluid 
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velocity,  effects  only  the  domain,  the  deflections  in  the  boundary  being 
only  a  secondary  effect  of  convection.  On  the  other  hand,  when  the 
upper  boundary  is  also  kept  at  a  fixed  temperature  we  have,  Ma  =  Bo  =  Cr 
=  0  and  so 

at  z  -  1,  w=Dw=e=o  (5.1-14) 

This  then  is  the  pure  Rayleigh-Benard  problem  with  all  the  nonlin- 
earities  only  in  the  domain  and  the  resulting  nonquiescent  solution  also 
manifests  itself  in  the  domain  as  convection.  The  Rayleigh-Marangoni 
problem  described  by  eqns.  (5.1-8)  -  (5.1-13)  is  a  mixed  problem  with 
nonli near i ties  in  the  domain  and  the  boundary  but  the  convective  solu- 
tion resulting  from  these  nonli near i ties  shows  up  mainly  in  the  domain. 

5.2.  The  Augmented  Morphological  Problem 

As  can  be  seen  from  Section  5.1,  in  the  Rayleigh-Marangoni  problem 
there  is  a  Ra-Ma  domain-boundary  duality  which  does  not  seem  to  exist  in 
morphological  instability.  From  the  problem  description  in  Chapter  3  it 
is  easy  to  see  that  all  the  nonlinearities  for  this  problem  lie  only  in 
the  liquid-solid  interface.  This  is  a  limitation  because  by  virtue  of 
being  on  the  boundary  the  Sekerka  number  is  unique  and  hence  also  has  a 
unique  eigenfunction  and  is  insufficient  when  solutions  to  inhomogeneous 
versions  of  the  linearized  morphological  instability  problem  are  needed, 
as  in  Chapter  6  where  "imperfections"  are  considered.  This  difficulty 
also  crops  up  in  the  pure  Marangoni  problem,  but  the  Rayleigh-Marangoni 
problem  comes  to  the  rescue,  as  there  are  countably  many  corresponding 
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values  of  Ra  for  each  value  of  Ma  and  hence  also  countably  many 
eigenf unctions  forming  a  complete  set  (see  Rosenblat,  Homsy  and  Davis 
(1981)).  The  naturally  occurring  duality  of  Ma  and  Ra  enables  solutions 
to  inhomogeneous  problems  to  be  obtained  in  a  straightforward  manner. 

In  morphological  instability  there  is  no  such  obvious,  naturally 
occurring  boundary-domain  duality  and  it  is  necessary  to  create  one.  To 
avoid  confusion  we  will  refer  to  the  pure  morphological  problem  of 
Chapter  3  as  the  Sekerka  problem,  and  (by  analogy  with  the  Rayleigh- 
Marangoni  problem)  set  up  an  eigenvalue  problem,  with  the  eigenvalue  in 
the  domain,  which  we  will  call  the  "augmented  morphological  problem." 


liquid:     V2p4  =  -  MPjl  (5.2-1) 

2       3q£ 
V  q£  "  V  JT  '   °  (5.2-2) 


solid:      V  pg  =  -  Mps  (5.2-3) 

2       9qs 
nV  qs  "  V  dT  '   °  (5.2-4) 


where  M  is  the  eigenvalue  which  we  label  as  the  morphological  number. 
The  boundary  conditions  are 

at  z  =  -1,     p4  =  q^  =  0  (5.2-5) 

at  z  =  s,      ps  =  0  (5.2-6) 

at  Z  -  0,      pa  -  ps  +  t(G£  -  G3)  =  0  (5.2-7) 
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p£  ♦  Seq^  *  tto     ♦  SeGlc)    +  A  1|  .  0 

3r 


(5.2-8) 


3P„     3P, 


9z 


9z 


(5.2-9) 


kq.  -  q  +  t(kG„   -  G  )  =  0 
X.    s      Ic         sc 


(5.2-10) 


3Qo    3q, 


3z 


n  5T  "  Vq*  +  vqs  =  ° 


(5.2-11) 


with  periodic  conditions  at  the  lateral  boundary  at  r  =  R.  We  can  now 
separate  variables  expressing  the  horizontal  dependence  as  zeroth  order 
Bessel's  functions  of  the  first  kind  JQ(air/R),  where  a,  are  the  zeros 
of  the  first  order  Bessel's  functions  of  the  first  kind  J,. 

00      00 

p.(r,z)  =  I   Z     p0..(z)  J  (a.r/R)       (5.2-12) 
*       i-1  j=i   "J     °  x 


In  addition  if  we  take  Fourier  transforms  in  the  horizontal  direc- 
tion and  solve  for  q. .  . ,  q  ..  and  t.  the  equations  reduce  to  a  system 

S.ij   sij      i  J 

in  p„ . .  and  p  . . .  If  we  define  a  column  vector 
Kiij     Ksij 


Q.. 


U 


'ilj 


'sij 


(5.2-13) 


and  a  matrix  differential  operator  L^ 


L. 

i 


(D* 


<•> 


BY.(D' 


>?> 


(5.2-14) 
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then  the  domain  equations  reduce  to 


L.Q. .  =  -  M. .Q. . 
i  ij     ij  ij 


(5.2-15) 


where  Y. 


(G.  +  Se  8 
I  o 


a. A  )/(G  +  Se  G     -  a2A  ) 
i      s       c    1  ' 


(5.2-16) 


Gc  "  (G*cB  +  Gsc)/(B  +  k) 


(5.2-17) 


B  = 


(a   -  v/2tanha„ . )tanha  .s 
x-i x,i si 

(na   +  v/2tanha  . s)tanha„. 
si  si      %i 


(5.2-18) 


au  -  Caf  .  v2/„"= 


(5.2-19) 


a3.    =  (af  ♦  vW)1/2 


(5.2-20) 


The  boundary  conditions  become 


ati--1,  p         =o 


(5.2-21) 


at   z   =  s,  p    .      =   0 

si  J 


(5.2-22) 


at  z  =  0, 


where 


B.Q.  . 

=   0 

B.    = 

1 

-Y, 

i 

D 

-6D 

(5.2-23) 


(5.2-24) 
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Defining  an  inner  product 

<Qij'Qik>  =  _|  Ply  Pii„  dz  +  I  Kij   ?sik  dz      (5'2"25) 

where  the  "*"  refers  to  the  adjoint  eigenf unction  and  "-"  the  complex 
conjugate. 

It  can  easily  be  seen  that  the  system  described  by  eqns.  (5.2-13)  - 
(5.2-24)  is  self-adjoint  in  this  inner  product  and  so  the  eigenfunctions 
Q^.:  are  complete.  Solving  the  system  we  get 


*ij 


A..  .Sin  A.  .-  a2  (z  +  1) 


A  . .Sin/M. .-  a2  (s  -  z) 
sij     ij   i 


(5.2-26) 


where  M« .   are  solutions  of  the  equation 


Se  =   (-  gt  +  a.A    )/  GQ 


(5.2-27) 


k.G.tan/M. .-a.    S  +  k  G  tan/M. .-  a. 

G-   -      Li  LI        * 3   S      .       1J    .   1  (5.2-28) 

k.tan/M. .-  a     S  +  k  G  tan/M. .-  a 
i  IJ        1  ss  IJ         1 


Here  A.. .    and  A   . .    can  be  determined  from  the  normalizing  condition 
fcij  sij 


<Q.  .,Q.,  >   =   5., 

ij      ik  jk 


(5.2-29) 


where  6. .  is  the  Kronecker  delta. 
IJ 
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5.3.  Comparison  of  Morphological  Instability  with 
Rayleigh-Marangoni  Convection 


The  principal  aim  of  this  section  is  to  relate  the  mathematical 
characteristics  of  the  two  problems  so  that  we  may  introduce  some  of  the 
extensive  mathematical  techniques  used  to  study  Rayleigh-Marangoni 
convection  to  morphological  instability,  but  we  will  make  some  physical 
comparisons  as  well.  The  augmented  morphological  problem  described  in 
Section  4.2  is  similar  to  Rayleigh-Marangoni  convection.  The  augmented 
problem  is  self  adjoint  but  the  Rayleigh-Marangoni  problem  is  nonself 
adjoint.  Both  have  an  infinity  of  eigenvalues  and  corresponding  eigen- 
f unctions,  but  while  completeness  of  the  eigenspace  is  assured  for  the 
former,  special  theorems  are  required  to  show  this  for  the  latter  (cf. 
Nadarajah  and  Narayanan  (1987)).  It  should  also  be  noted  that  while  the 
Rayleigh-Marangoni  problem  attempts  to  describe  a  realistic  situation, 
the  augmented  morphological  problem  was  artificially  created  in  order  to 
solve  inhomogeneous  versions  of  the  Sekerka  problem  described  in  Sec- 
tions 3.1  and  3.2. 

This  brings  us  to  the  question  whether  there  is  a  practical  situa- 
tion which  is  described  by  this  mathematical  concoction.  The  difficulty 
in  coming  up  with  one  stems  from  another  important  difference  between 
the  two  problems.  In  the  pure  Rayleigh-Benard  problem  (where  Ma,  Bo  and 
Cr  are  all  zero)  the  nonlinearity  is  in  the  domain  and  the  instability 
too  manifests  in  the  domain  as  convection.  Even  in  the  pure  Marangoni 
problem  (where  Ra  is  zero)  where  the  nonlinearity  is  in  the  boundary, 
the  instability  is  still  mainly  in  the  domain.  In  contrast,  in  the 
Sekerka  problem  the  nonlinearities  and  the  resulting  instability  show  up 
in  the  boundary.   Though  there  are  other  boundary  effects  (like  kinetic 
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undercooling)  which  can  cause  morphological  instability  the  only  domain 
effect  which  could  give  M  physical  significance  is  a  heat  source  term  in 
the  form  MT  or  Me~E/RT  (see  Joseph  (1965)).  We  do  not  know  of  any 
experiment  where  morphological  instability  was  observed  as  a  result  of  a 
heat  source  in  the  melt  or  the  crystal  but  if  one  does  exist  it  will 
provide  the  true  analogy  to  Rayleigh-Benard  convection. 

This  is  relevant  as  Hurle  (1985)  has  attempted  a  comparison  between 
the  Rayleigh-Benard  problem  and  the  Sekerka  problem.  It  can  now  be  seen 
that  the  Sekerka  problem  can  only  be  compared  to  the  pure  Marangoni 
problem,  with  Se  corresponding  to  Ma  and  A  corresponding  to  the  reci- 
procal of  Bo.  Besides,  for  periodic  lateral  boundary  conditions,  the 
eigenvalues  of  both  problems,  Se  and  Ma,  are  unique. 

Based  on  this  comparison  we  can  make  an  important  conjecture. 
Vrentas,  Narayanan  and  Agrawal  (1981)  have  shown  that  for  the  Marangoni 
and  the  Rayleigh-Marangoni  problems  when  the  nonperiodic  no-slip  condi- 
tion for  velocity  is  imposed  at  the  sidewalls,  the  eigenvalue  Ma  is  no 
longer  unique  and  has  countably  many  values.  In  other  words  when  the 
walls  are  a  finite  distance  apart  Ma  has  many  values,  but  as  they  are 
gradually  moved  apart  we  have  "spectral  crowding"  and  in  the  limit  when 
they  are  sufficiently  far  apart  to  impose  the  periodic  boundary  condi- 
tion of  total  slip,  all  the  values  of  Ma  coalesce  into  a  unique  num- 
ber. Recently,  following  Coriell  et  al.  (1980),  several  workers  have 
looked  at  the  coupled  problem  of  morphological  instability  with  solutal 
convection  and  all  have  assumed  periodic  boundary  conditions.  We  sus- 
pect that  here  too  if  the  no-slip  condition  for  velocity  at  the  side 
wall  is  imposed,  the  Sekerka  number  will  no  longer  be  unique.  All  this 
raises  the  question  of  completeness  of  the  Marangoni  and  the  Sekerka 
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eigenspaces  and  its  probable  that  generalized  eigensolutions  (see  Nai- 
mark  (1967))  are  needed  when  Ma  and  Se  are  chosen  as  eigenvalues. 

When  these  two  problems  are  considered  in  a  finite  container  we  can 
see  yet  another  difference.  Both  problems  have  simple  eigenvalues 
except  at  certain  aspect  ratios  of  the  container  where  two  horizontal 
modes  can  coexist  (cf.  Rosenblat,  Homsy  and  Davis  (1982)).  In  a  typical 
experiment  of  Rayleigh-Benard  convection  we  would  expect  to  see  a  dozen 
or  so  convection  cells  (see  for  example  Koschmieder  (1967))  and  increas- 
ing or  decreasing  the  number  of  cells  by  one  can  significantly  effect 
the  problem.  Hence  the  multiple  points  in  this  problem  are  extremely 
important  and  have  been  the  subject  of  study.  But  in  morphological 
instability  a  single  alloy  crystal  can  contain  hundreds  of  individual 
cells  and  the  addition  or  loss  of  one  has  hardly  a  noticeable  effect  on 
the  problem  and  consequently  multiplicity  of  the  lateral  eigenfunctions 
loses  its  significance.  In  addition  unlike  the  Rayleigh  number  it  is 
well  known  that  near  the  critical  wave  number  aN,  the  critical  value  of 
the  Sekerka  number  SeQ  hardly  changes  (see  for  example  Coriell,  McFadden 
and  Sekerka  (1985))  and  the  choice  of  aN  has  very  little  effect  on 
SeQN.  Conversely,  the  choice  of  the  operating  Se  will  have  a  tremendous 
impact  on  the  resulting  wavenumber  (cf.  Ramprasad  and  Brown  (1987)). 

Other  differences  have  been  mentioned  in  Chapter  4.  Both  the 
Marangoni  problem  and  the  Sekerka  problem  have  symmetric  bifurcation 
diagrams  near  the  bifurcation  point  for  two-dimensional  rolls  and  rec- 
tangular cells  and  nonsymmetric  curves  for  hexagonal  cells  and  cylindri- 
cal rolls.  But  in  morphological  instability  the  curves  can  be  "forward 
bending"  or  "backward  bending"  depending  on  the  operating  conditions, 
whereas  in  Marangoni  convection  the  curves  are  forward  bending  every- 
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where.   Hence  the  occurrence  of  subcritical  instability  is  more  wide- 
spread in  morphological  instability. 


CHAPTER  6 
BIFURCATION  BREAKING  IMPERFECTIONS 


6.1  .  Nature  of  Imperfections 

When  the  morphological  instability  problem  was  formulated  in  Chap- 
ter 3,  several  effects  were  ignored  and  the  resulting  problem  is  an 
idealized  or  "perfect"  one.  Inclusion  of  these  can  alter  the  problem  in 
several  ways,  for  example  kinetic  undercooling  of  the  melt  becomes  an 
important  effect  in  rapid  solidification,  but  all  it  does  is  alter  the 
onset  condition  for  morphological  instability.  In  the  parlance  of 
bifurcation  theory,  an  "imperfection"  is  an  effect  on  the  "perfect" 
problem  which  alters  it  in  a  specific  way.  Such  an  imperfection  will 
cause  the  morphological  instability  problem  not  to  have  a  planar  solu- 
tion at  all  even  below  the  onset  condition,  and  these  are  known  as 
bifurcation  breaking  imperfections. 

The  effect  of  a  typical  imperfection  on  the  bifurcation  diagram  is 
shown  in  Fig.  6-1  .  The  broken  line  is  the  solution  in  the  presence  of 
imperfection  and  it  can  be  seen  that  the  interface  will  be  nonplanar  for 
all  nonzero  values  of  Se.  Obtaining  solutions  to  the  problem  with 
imperfections  is  extremely  difficult  and  we  will  only  seek  asymptotic 
solutions.  Hence  the  problems  to  be  considered  should  have  very  small 
imperfections.  Under  such  conditions  the  method  of  matched  asymptotic 
expansions  of  Matkowsky  and  Reiss  (1977)  can  be  employed  and  here  it 
will  be  used  in  a  way  similar  to  the  work  of  Tavantzis,  Reiss  and  Mat- 
kowsky (1978)  for  the  Rayleigh-Benard  problem. 
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Se  =Seo  n 


Figure  6-1.   Imperfect  bifurcation  diagram  showing  inner 
and  outer  expansions 
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The  method. is  fairly  straightforward.  The  variables  are  expanded 
asymptotically  with  the  imperfection  parameter  about  the  planar  and  the 
nonplanar  solutions  and  two  outer  expansions  0Q  and  0,  are  obtained  as 
shown  in  Fig.  6-1 .  At  the  bifurcation  point  SeQN  these  expansions  break 
down  and  it  is  necessary  to  have  inner  expansions  I1  and  Ip  near  SeoN 
and  to  join  the  corresponding  0q  and  0^  matching  conditions  have  to  be 
specified. 

The  imperfections  that  can  be  analyzed  in  this  fashion  must,  of 
course,  be  small  effects  else  a  full-blown  nonlinear  solution  will  be 
needed.  Two  of  the  most  important  effects  which  are  habitually  ignored, 
namely  imperfect  insulation  of  the  ampoule  wall  and  advection  in  the 
melt,  are  such  imperfections  and  readily  lend  themselves  to  this  type  of 
analysis.  In  Chapter  3  when  the  morphological  instability  problem  was 
modelled  by  a  uniform  formulation,  it  was  mentioned  that  one  reason  for 
this  was  to  consider  a  finite  crystal/melt  region.  This  finiteness  was 
only  in  the  vertical  direction  and  in  the  horizontal  direction  the 
imposition  of  periodic  boundary  conditions  effectively  meant  that  the 
ampoule  side  walls  were  infinitely  far  apart.  The  two  imperfections 
that  are  to  be  considered  are  caused  by  nonperiodic  lateral  boundary 
conditions  and  another  way  of  looking  at  the  effect  of  these 
imperfections  is  to  say  that  the  container  is  now  being  considered  to  be 
finite  in  the  lateral  as  well  as  the  vertical  direction. 

6.2.   Imperfection  Due  to  Heat  Loss 

As  mentioned  in  the  last  section,  it  has  been  customary  in  this 
problem  to  assume  periodic  boundary  conditions  laterally,  which  is 
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equivalent  to  assuming  that  the  walls  of  the  ampoule  are  perfectly 
insulated  or  that  they  are  so  far  apart  that  their  effect  can  be  ig- 
nored. In  practice  neither  of  these  is  likely  to  be  achieved  and  here 
we  will  examine  the  effect  of  a  small  amount  of  heat  loss  or  heat  gain 
from  the  wall  on  the  problem.  If  we  take  the  ampoule  to  be  cylindrical 
with  radius  R, 


3T 
r  =  R,  liquid:     -k£  |~  -  6f£(z)  (6.2-1) 

dT 
SOlid:      ~ks  3T"  =   6fs(z)  (6.2-2) 


where  t%   and  fg  are  such  that  f %{0   =  f  (5)  and  f£(-1 )  =  f  (s)  =  0. 
If  f^  and  fg  are  positive  it  will  mean  heat  loss  and  if  they  are  nega- 
tive, heat  gain.   If  we  make  the  transformations 


f.Cz) 
T£(z,r)  =  5  — r  +  e^z.r)  (6.2-3) 

f  (z) 
and  T  (z,r)  =  <5  — r  +  0  (z,r)  (6.2-4) 

3  K  S 


and  substitute  these  in  the  steady-state  versions  of  eqns.  (3.1-15) 
(3.1-24),  the  temperature  equations  in  the  domain  become 


V2e£  =  -  5rD2f£/k£  (6.2-5) 


V26s  =  -  5rD2f  /k  (6.2-6) 
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The  outer  boundary  conditions  will  remain  unchanged  but  at  the  interface 
eqn.  (3.1-21)  becomes 

\  ~  9S  ■ "  6r(Vk* "  W  (6-2"7) 


9^  «■  SeC£  +  A  H     =  -  -Srf^/k^         (6.2-8) 


and  eqn.  (3.1-22)  converts  to 


V6  -n  -  BV6  -n  =  Lv  -  6r(Df  -  Df  )        (6.2-9) 

«*  o  X/      3 


In  order  to  solve  this  system  we  will  be  treating  the  heat  loss  as 
an  imperfection  on  the  perfect  problem  (i.e.  when  6=0).  The  perfect 
problem  is  of  course  the  Sekerka  problem  of  Section  3.1. 

6.3.   The  Outer  Expansions 

As  this  problem  has  been  defined  in  a  finite  geometry,  the  number 
of  cells  are  fixed  by  the  container  size  and  the  growth  conditions.  We 
will  choose  these  so  that  aN  is  very  close  to  amin  the  wavelength  corre- 
sponding to  SeQmin  the  least  value  of  SeQN.  Another  important  decision 
is  the  selection  of  the  wave  pattern  and  our  analysis  is  done  for  cylin- 
drical rolls.  An  objection  to  this  could  be  raised  on  the  grounds  that 
in  most  experiments  it  is  the  hexagonal  pattern  which  is  observed.  We 
justify  our  assumption  by  once  again  making  a  comparison  with  Rayleigh- 
Benard  convection.  It  appears  that  in  Rayleigh-Benard  convection  the 
wave  pattern  selection  is  strongly  influenced  by  the  container  size  and 
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shape,  with  the  hexagonal  pattern  prevailing  for  all  shapes  in  wide 
containers  while  in  narrower  ones  the  container  shape  determining  the 
pattern  (e.g.,  cylindrical  rolls  for  circular  containers.  See  Kosch- 
mieder  (1967)).  Since  the  principal  aim  of  this  paper  is  to  study  a 
finite  geometry  effect  on  morphological  instability,  the  cylindrical 
roll  pattern  would  be  the  logical  choice.  This  is  especially  valid  for 
experiments  such  as  those  of  Peteves  (1986)  where  ampoules  of  radius 
0.025  inches  were  used. 

In  Chapter  4  it  was  shown  that  the  bifurcation  diagram  is  unsymme- 
tric  for  cylindrical  rolls  and  the  form  of  the  outer  expansions  0Q  and 
0.]  are  shown  schematically  in  Fig.  6-1.  If  we  use  superscript  o  to 
identify  the  problem  with  perfect  insulation  we  can  seek  solutions  by 
means  of  asymptotic  expansions  about  the  perfect  problem. 


L  =  Z     ej  6k  (6.3-D 

*   k=0  * 


with  similar  expansions  for  6  ,  C. ,  C  and  c.   Substituting  these  expan- 

5     X*     5 

sions  into  the  problem  and  collecting  the  terms  of  order  6  we  get  an 
inhomogeneous  linear  system.  If  we  separate  variables  horizontally 


e!  -  I     el  J  (a.r/R)  (6.3-2) 

%        .    ,      % .    o  1 
i  =  1   i 

and  eliminate  C..  ,  C  and  z.     we  get  for  the  expansion  about  the  planar 
solution 


L.*1  =  f1  (6.3-3) 

li    i 
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where 


*.    = 

1 


31 

s. 


(6.3-4) 


t\  -  -  v2 


Vk* 


f  /k 

3       S 


(6.3-5) 


:.    =  J    r2J   (a.r/R)dr/   J    rJ2(a. 
1        ;  oi  *        o     1 


r/R)dr 


(6.3-6) 


The  boundary  conditions  are 


at  z  =  -1 , 


41 


=  0 


(6.3-7) 


at  z  =  s , 


si 


=  0 


(6.3-8) 


at  z  =  0, 


11    l 


(6.3-9) 


where 


I. 

i 


f./k.     -Y.f  /k 
%     %  l  s  s 


Df, 


-Df 


(6.3-10) 


The  eigenvalue  problem  of  Section  5.2  has  been  shown  to  have  a 
complete  set  of  eigenfunctions  and  hence  can  be  used  to  solve  the  above 
system. 
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00      00 


♦  =  I        I     <0  ,  *J>  Q  J  (a.r/R)        (6.3-11) 
1-1  j=1    J    x        1J  °  1 


1     j(Qii'*xi) 

<Q.  .,#'>  = sJ — L. 

ij   l  M.  . 


<Q. .,f]> 

z=0      ij 


(6.3-12) 


where  J(Q_,h.  )L=0  is  the  bilinear  concomitant  evaluated  at  z=0 


J(W|Z«0  "  "  P*i(Df*  "  DVVk&  +  DPU(fH/k£  "  VA}Ii 

(6.3-13) 

If  we  assume  that  h^  and  f  are  nonzero,  then  when  Se  =  Se  N,  MN1 
becomes  zero  and  the  outer  expansion  (6.3-11)  -  (6.3-13)  will  fail. 

When  we  expand  about  the  nonplanar  solution  we  end  up  with  a  much 
more  complicated  system.  But  since  we  have  expressed  the  solution  to 
the  perfect  nonplanar  problem  itself  in  terms  of  a  perturbation  series 
in  e,  we  can  do  the  same  for  the  imperfect  problem. 


1    1     1 
*  =  *Q  +  £*1  +  ...  (6.3-14) 


If  we  take  the  zeroth  order  problem  and  once  again  separate  vari- 
ables laterally, 

CO 

*o  =  .E  *oi  J0(air/R)  (6.3-15) 

and  eliminate  C.  . ,  C   .  and  c  .  ,  we  get  the  same  system  described  by 

X.01     SOI        01 

eqns.  (6.3-3)  _  (6.3-10)  but  with  #  .  as  the  variable.   Again  using  the 
eigenfunctions  of  Section  5.2,  the  solution  is 


82 


o»  eo 


•o'i!,  ji,  <Qia'*oi>Qij  w/R> 


(6.3-16) 


with  <Q. j, #Qi>  taking  the  same  form  as  eqns .  (6.3-12)  and  (6. 3-13). 
Since  these  equations  are  valid  only  for  Se  close  to  SeQN,  we  can  write 


an  expansion  for  M 


N1 


MN1(Se)  =MN1(SeoN)  ♦  (Se  -  SenM) 


dM 


N1 


oN  dSe 


Se  =  Se 


ON 


(6.3-17) 


which  simplifies  to 


MN1(Se)  =  0  ♦  eSe^  ♦  o<«') 


(6.3-18) 


where 


dM 


M'   = 


N1 


N1    dSe 


Se  =  Se 


oN 


(6.3-19) 


and  Se1N  has  been  derived  in  Section  4.2.   Hence  the  outer  expansion 
about  the  nonplanar  solution  is 


*  =  *c   +    ~*oN  Jo(aNr/R)    +   0(e    >    + 


(<Q^    ,fj>   -   J) 

5[        £Se1NM-N1        QN1    Jo(aNr/R)    +   0(eO)]    +   0(6    } 


(6.3-20) 
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6.4.  The  Inner  Expansions 

The  outer  expansions  described  in  the  previous  section  fail  at  Se  = 
SeoN'  squiring  inner  expansions  in  this  vicinity.  For  Se  very  near 
SeoN 


CO 

Se(u)  =  Se   +  5yb  ♦  I     5kukb/k!       (6.4-1) 

k=2 


«(u)  =  yc  (6.4-2) 


where  b  and  c  are  integers  to  be  determined  such  that  the  solution  does 
not  fail  at  Se  =  SeQN.  We  can  now  expand  the  solutions  in  orders  of  p. 


u(r,z,u)  =  9(r,z,Se,6)  =  I     uk(r,z)yk/k!  (6.4-3) 

k=0 


v(r,z,u)  =  C(r,z,Se,6)  =  Z     v  (r,z)y  /k!  (6.4-4) 

k=0 

CO 

w(r,z,y)  =  s(r,Se,5)  =  Z  wk(r)uk/k!  (6.4-5) 

k=0 


Here  the  superscripts  on  u,  v  and  w  are  not  powers  but  indices. 

Substituting  these  expansions  into  the  problem  set  up  in  Section  6.2  and 

collecting  the  terms  of  order  \i,   we  once  again  get  an  inhomogeneous 

problem.   Separating  variables  and  eliminating  v  ,  v  and  w,  the  problem 

x<   s 

reduces  to 


L.*1  =  6'f!  (6.4-6) 

ill 


with  boundary  conditions 
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at  z  =  -1 , 


v  -  ° 


(6.4-7) 


at  z  =  s 


u  .  =  0 

31 


(6.4-8) 


at  z  =  0, 


BA  -  5'hJ 


(6.4-9) 


where 


(6.4-10) 


and 


.,   d6 

6  =dH 


u=0 


Se' 


dSe 
dy 


»«0 


(6.4-11) 


The  variables  Lt   and  Bi  have  been  defined  in  Section  5.2  and  f1 

1  1 

and  hi   in  Section  6.3.   But  Y,  now  becomes 


Yi  =  (G£  +  SeoN  Gc  "  aiA)/(Gs  +  SeoN  Gc  '  4A)  ^^'^ 


This  means  that  the  operator  Lj  is  singular  when  i=N,  j=1  and  hence 


,J\ 


for  solutions  to  exist  6  f N  has  to  orthogonal  to  the  null  space  of  L^ 
Using  this  solvability  condition  with  QN1  we  get 


6'[J(QNT^ 


z=0 


<QN1'V3  =  ° 


(6.4-13) 
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From  this  we  can  conclude  that  6'  =  0  (assuming  J  *  <Q  fl>)  and 
that  c>1.  This  makes  the  system  (6.4-6)  -  (6.4-12)  homogeneous  and 
so  4>  are  constant  multiples  of  *  . 


*1  =  A* 


(6.4-14) 


Proceeding  to  the  next  order  we  have 


2       1 
L.*  =  6"f 
11       l 


(6.4-15) 


and  at  z=0, 


2    2 
B.*7  =  h 

li    l 


(6.4-16) 


h.  -  2ASe' 
l 


C  .  G  /(G.  +  Se  „  G  -a7A  ) 
oi  c  I  oN  c  i 


+  2A  I.,Q 

i1 


(6.4-17) 


where 


[.,  =  J  rJ3(a. 
i1   i   o  i 


r/R)dr/  J  rJ  (a.r/R)dr 
0   °  x 


(6.4-18) 


and  Q  is  a  complex  function  of  the  linear  stability  variables.   The 
solvability  condition  for  i=N,  j=1  now  gives 


J(QN1,h^)  -6"   <Qm,tl>   =  0 


(6.4-19) 


Hence  we  can  take  &' '  *  0  and  Se'  *  0  and  so  b=1  and  c=2.  Equation 
(6.4-19)  is  a  quadratic  in  A  with  two  real  roots  A1  and  A2  of  opposite 
signs,  corresponding  to  two  inner  expansions  I-|  and  ^^ 
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Se  =  SeQN  +  £S1/2  +  0(6)  (6.4-20) 


#  -  f  °  +  Av*q51/2  +  0(6)  (6.4-21) 


In  terms  of  y,   the  outer  expansion  0Q  about  the  planar  state  of 
eqn.  (6.3-11)  can  now  be  expressed  by 


#il  "  ^NT^  "  J)QN1  Jo(aNr/R)/5MN1      (6-4"22) 


and  the  nonplanar  expansion  0,  of  eqn.  (6.3-20) 


#il  =  ^^1'^  "  J)/CMN1  +  DN^/SelN]QN1  Jo(V/R)     (6'i|-23) 


where  we  expressed  *qN  as  DNQN1  J  (a  r/R) .   In  attempting  to  match  these 
with  the  inner  expansions  of  eqn.  (6.4-21)  we  get 


lim  A1  =  lim  A2  -«Q^  ,fj>  -  J)/5M'       (6.4-24) 


lim  A1  =  lim  A2  =  D  5/Se  (6.4-25) 

£-*■  00  £->   —00 


Hence  using  A1  in  eqn.  (6.4-21)  will  give  the  inner  expansion  I1 

for  matching  0Q  for  Se  <  SeQN  with  01  for  Se  >  SeoN.  A2  will  result  in 

I2,  matching  0Q  for  Se  >  SeoN  with  01  for  Se  <  SeQN  as  shown  in  Fig.  6- 

1. 
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6.5.  Imperfection  Due  to  Advection  in  the  Melt 

In  modelling  morphological  instability  during  solidification  it  is 
customary  to  neglect  the  difference  in  density  between  liquid  and 
solid.  Though  this  difference  is  very  small  (e.g.  for  the  Pb-Sn  system, 
the  ratio  of  densities  ps/p£  is  1.035)  the  volume  contraction  upon 
solidification  results  in  closed  streamlined  flow  in  the  melt  which 
fundamentally  alters  the  planar  state  of  this  problem.  Obtaining  the 
solution  to  the  planar  state  with  this  flow  is  difficult  enough  even 
without  attempting  the  formidable  task  of  solving  the  nonplanar  prob- 
lem. In  order  to  include  the  effect  of  the  flow  on  morphological  insta- 
bility previous  workers  have  therefore  relied  on  simplifying  assumptions 
or  looked  at  limiting  cases. 

The  scenario  in  a  typical  solidification  experiment  is  schemati- 
cally shown  in  Fig.  6-2.  The  heating  coils  maintain  constant  tempera- 
tures at  positions  z  =  -I  and  z  =  s  in  the  ampoule.  As  solidification 
proceeds,  the  ampoule  is  pulled  in  the  positive  z-direction  at  a  velo- 
city v  in  order  that  the  liquid-solid  interface  will  remain  stationary 
at  z  =  0.  To  fill  the  space  created  by  volume  contraction  upon  solidi- 
fication, the  melt  will  move  towards  the  crystal  with  a  bulk  velocity 
of  v(ps/p^  -  1)  and  this  process  is  commonly  referred  to  as  "advec- 
tion." In  the  presence  of  the  rigid  ampoule  walls,  this  flow  also 
resembles  the  bolus  flow  of  slugs  through  a  pipe,  resulting  in  closed 
streamlined  flow  in  the  melt. 

In  the  literature  we  find  three  approaches  to  tackling  this  prob- 
lem, the  least  of  which  is  the  work  of  Caroli,  Caroli,  Misbah  and  Roulet 
(1985b)  who  ignored  the  rigid  side  walls.   Then  the  only  effect  of 
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Figure  6-2.   Crystal  growth  with  advection  in  the  melt 
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advection  is  a  negligible  modification  of  the  growth  velocity  in  the 
melt.  The  other  two  approaches  are  more  substantial  and  analize 
limiting  cases  of  this  phenomenon.  Since  the  traditional  formulation  of 
morphological  instability  concentrates  only  at  a  "boundary  layer"  region 
of  the  liquid-solid  interface,  in  such  a  model  the  closed  streamlined 
flow  shown  in  Fig.  6-2  can  be  approximated  by  Couette  flow  in  Region  I 
and  stagnation  flow  in  Region  II.  Delves  (1968,  1971  and  1974)  and 
Coriell,  McFadden,  Boisvert  and  Sekerka  (1984)  looked  at  the  effect  of  a 
forced  plane  Couette  flow  and  showed  that  the  onset  of  morphological 
instability  is  somewhat  suppressed  for  disturbances  in  the  flow 
direction,  while  disturbances  perpendicular  to  the  flow  were 
uneffected.  Recently  McFadden,  Coriell  and  Alexander  (1988)  examined 
the  effect  of  a  plane  stagnation  flow  on  disturbances  perpendicular  to 
the  flow  and  here  too  the  flow  was  found  to  be  stabilizing. 

Based  on  these  two  limiting  cases  it  might  be  tempting  to  conclude 
that  advection  generally  stabilizes  the  liquid-solid  interface.  In  this 
and  the  next  few  sections  we  will  embark  upon  a  more  complete  analysis 
of  advection  than  has  been  attempted  so  afar  and  show  that  the  above 
assertion  is  questionable  at  best.  In  our  model  we  consider  morphologi- 
cal instability  with  advection  in  the  absence  of  natural  convection.  We 
choose  not  to  include  natural  convection  as  we  wish  to  study  the  effect 
of  advection  only  and  natural  convection  will  only  further  complicate  an 
already  nontrivial  problem.  Besides  it  has  been  shown  very  conclusively 
by  Coriell,  Cordes,  Boettinger  and  Sekerka  (1980)  and  by  Caroli,  Caroli, 
Misbah  and  Roulet  (1985a)  that  the  convective  and  morphological  modes 
are  decoupled  except  at  points  where  they  become  comparable.  Hence  our 
model  will  be  valid  for  the  high  growth  velocity  region  where  morpholo- 
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gical  instability  is  dominant  and  also  in  low  gravity  environments  where 
natural  convection  can  be  neglected. 

The  experimental  set  up  was  briefly  described  earlier.   The  steady 
state  domain  equations  in  dimensionless  form  are 


V2T£  =  0  (6.5-1) 


A.      3C* 


V  Ci  ~   V  W  ~  U'VCi   -   °  (6'5"2> 


V2Ts  =  0  (6.5-3) 


2       8C 
nV  Cs  "  V  W~   =  °  (6.5-4) 


Here  u  is  the  velocity  of  the  closed  streamlined  flow  in  the  melt 

caused  by  advection  and  5  is  (Pg/P^  -  1).  The  boundary  conditions  are 

at  z  =  -1,     Tl   =  T1 ,  C^  =  C1  (6.5-5) 

at  z  =  s,      Ts  =  T2  (6.5-6) 

at  z  -  5,      T£  =  Ts  =  -  SeC^  +  A  H  (6.5-7) 

VT  -n  -  3VT  «n  =  Lv  (6.5-8) 
x»       s 


kC&  =  Cs  (6.5-9) 


VC  -n  -  (1+5)vC.  =  nVC  -n  -  vC  (6.5-10) 

X»  x,  s        s 
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Before  we  can  write  the  equations  for  u,  further  assumptions  are 
necessary.  The  domain  of  the  velocity  equations  extend  over  the  entire 
melt  region,  not  just  the  depth  I  that  we  have  considered.  In  the  float 
zone  technique  this  melt  region  will  still  have  a  constant  depth  as 
solidification  proceeds  but  in  Bridgman  growth  the  melt  region  will 
continuously  shrink.  Since  crystal  growth  velocities  are  usually  so 
small,  even  in  the  latter  technique  it  takes  awhile  before  there  is  an 
appreciable  change  in  the  melt  depth.  Hence  over  a  short  time  span  a 
constant  depth  assumption  will  be  valid  even  for  Bridgman  growth.  We 
will  also  assume  that  the  melt  is  bounded  on  all  sides  by  rigid  walls. 
This  will  require  a  container  for  zone  refining  and  a  very  viscous 
encapsulant  for  the  melt  in  Bridgman  growth. 

Under  these  assumptions  the  velocity  problem  has  been  solved  by 
Duda  and  Vrentas  (1971)  for  low  velocities  and  here  we  will  only  give 
the  solution  and  refer  the  interested  reader  to  their  paper  for  the 
details . 


u  =  (5vR/r)3i|i/3zf    u  =  6v(1  -  (R/r)3i|>/3r)  (6.5-11) 


i|i  =  I     A  F  (r)Sinmr(l  +  z/Y)  +  I     G  (z)rJ,(b  r/R)     (6.5-12) 
n  n  n     1  n 

n=1  n=1 

F   (r)    =   [rl   (nirr/Y)I   (nir/h)    -  r2I]  (mr/hJIgdnrr/YVRj/IgCnir/h)      (6.5-13) 


G  (z)    =  2B   [Sinhb   (z+Y)/R  -    (z/Y+1 )exp(-b  z/R)Sinhb  h] 

n  nL  n  *       n  n   J 


+  2C   (z+Y)/R  expb  h  Sinhb  z/R  (6.5-14) 

n  n  n 
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p  O 

gn  =  J0(bn)f4bnhexpbnh  +  (2exPbnh  "  exP("bnh)  "  exP3t>nh)/h]         (6.5-15) 


Bngn  =  -HbnhexPbnh     I     A  Q        ♦  2(1    -  exp2b  h)      I     A  Q      (-1 )m        (6.5-16) 

m=1  m=1  m   um 


nS  C  =  Mb  hexpb  h  -  Sinhb  h)  I  A  Q 
n  n     n   *  n        n  '   ,  m  mn 

m=1 


(4b  h+2  -  2exp2b  h)   Z  A  Q   (-1  ) 
n  n       m  mn 

m=1 


ra 


(6.5-17) 


bn  =  £  nir[I0(nir/h)I2(nir/h)   _  lf(nw/h)]/l2(nir/h)  (6.5-18) 


2bn(mir/h)2I^(mir/h)J2(b  ) 
Qmn  =  ~~2 "T 2 2T2  (6.5-19) 


Y  =  im/l,  h  =  Y/R  (6.5-20) 


where  Jl^  is  the  total  melt  depth,  R  the  dimensionless  radius  of  the 
ampoule,  ur  and  uz  the  radial  and  axial  components  of  the  velocity,  J 
the  Bessel  function  of  the  first  kind  of  order  p  and  I  the  modified 
Bessel  function  of  the  first  kind  of  order  p.  Duda  and  Vrentas  have 
determined  the  first  twenty  coefficients  of  A  for  various  values  of  h 
and  we  will  use  these  in  our  calculations. 
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6.6  Nonexistence  of  the  Planar  State 

The  key  assumption  made  in  the  two  previous  approaches  to  solving 
the  advection  problem  was  that  a  stable  planar  state  was  still  possible 
even  in  the  presence  of  melt  flows.  But  this  is  true  only  if  the  velo- 
city profile  very  close  to  the  interface  is  assumed  to  be  in  a  special 
form,  the  vertical  component  of  velocity  being  independent  of  any  hori- 
zontal coordinates  and  the  horizontal  component  of  velocity  at  most  a 
linear  function  of  that  horizontal  coordinate.  As  can  be  seen  by  eqns. 
(6.5-11)  -  (6.5-19)  these  assumptions  can  only  be  viewed  as  limiting 
cases  and  in  general  ur  and  uz  will  be  complex  functions  of  both  r  and 
z.  In  this  section  we  will  use  this  to  prove  the  nonexistence  of  planar 
solutions  for  this  problem.  This  proof  can  also  be  used  to  show  the 
nonexistence  of  planar  solutions  for  the  case  of  imperfect  insulation 
considered  in  Sections  6.6  -  6.4. 

In  Section  3.1  where  we  assumed  that  6  was  zero,  a  planar  solution 
was  possible  with  the  temperatures  and  concentrations  functions  of  z 
only.  Then  the  energy  balance  condition  at  the  interface  (given  by  eqn. 
(6.5-8))  could  be  used  to  determine  the  constant  growth  velocity  v.  In 
other  words  the  choice  of  T1 ,  T2  and  C-,  fixed  v.  For  the  case  of  non- 
zero 6  we  will  show  that  planar  solutions  do  not  exist  by  contradic- 
tion. If  there  is  a  possible  planar  solution  eqn.  (6.5-2)  specifies 
that  C^  cannot  have  radially  independent  solutions.  The  solution  will 
be  in  the  form 


C.(r,z)  =  I     C0    (z)  J  (a  r/R)  (6.6-1) 

X,  sen     o  n 

n=1 


n 
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where  JQ  is  the  zeroth  order  Bessel  function  of  the  first  kind  and  a 
are  zeros  of  the  first  order  Bessel  function  of  the  first  kind.  Then 
the  boundary  conditions  (6.5-7),  (6.5-9)  and  (6.5-10)  will  require 
that  T^,  Ts  and  Cs  be  expressed  in  similar  expansions  in  terms  of  Bessel 
functions.  If  we  now  return  once  more  to  eqn.  (6.5-8)  to  calculate  v, 
we  can  see  that  v  too  should  be  expressed  in  terms  of  Bessel  functions 
radially,  but  this  contradicts  our  assumption  of  a  planar  solution. 
Hence  in  the  presence  of  advection  there  can  be  no  planar  solutions  to 
this  problem. 

6.7.  Asymptotic  Solution 


A  complete  nonplanar  solution  of  the  equations  described  in  Section 
6.6  will  require  tedious  numerical  calculations.  But  such  calculations 
may  be  unnecessary,  for  as  mentioned  before  for  most  systems  6  is  an 
extremely  small  number  and  hence  an  asymptotic  solution  should  be  pos- 
sible about  the  6=0  solution.  This  would  be  similar  to  the  method 
used  for  analyzing  heat  loss  at  the  ampoule  wall  and  here  too  we  will 
show  that  advection  is  a  bifurcation  breaking  imperfection  on  the  mor- 
phological instability  problem. 

But  unlike  the  problem  with  imperfect  insulation,  here  the  addi- 
tional terms  are  in  the  concentration  equations,  not  in  the  temperature 
ones.  Hence  as  done  earlier  we  cannot  eliminate  the  concentrations  and 
reduce  the  problem  to  only  the  temperature  equations  and  so  the  eigen- 
value problem  defined  in  Section  5.2  has  to  be  modified  to  include  the 
concentrations  in  the  column  vector.  Correspondingly  the  differential 
operator  Lj,  the  boundary  operator  B^^  and  the  inner  product  will  have  to 
be  changed  to  those  used  in  Chapter  3. 
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L.Q. .  =  -  M. .Q. . 

i  1J      ij  ij 


(6.7-1) 


L.  = 


2   2 
(D  -ap 


0 
0 


S*  <»2-v0-af) 


0 
0 


(D2-af) 


f2  <nD2-vD-naf) 


(6.7-2) 


<•*.♦>  -  ^  (TjT£  ♦  cjc£)dz  ♦  J  CT*Tg  ♦  Cs\)dz 


(6.7-3) 


B 


(1-Y.) 

D 

x/Y. 

l 

-vx/Y. 


Se 
0 


(xSe/Y.-k) 


[D-v(xSe/Y.+1-k)] 


L 

■3D 

0 
0 


0 
0 

1 

■nD 


(6.7-4) 


where  *  has  been  expanded  to  (T  ,  C. ,  T  ,  C  )  and 

JC     X»     3     S 


Yi  ■  {Gi  +  seGac  -arA)/(vGs) 


(6.7-5) 


x  =  (kG.   -  G  )/(G0  -  G  ) 
lo         sc    I         s 


(6.7-6) 


The  eigenvalues  Mj i  are  the  roots  of 


Se  =  (-  gt  +  a.A  )/  Gq 


(6.7-7) 
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GT  is  the  same  as  that  defined  in  eqn.  (5.2-8)  but  Gc  becomes 


Gc  =  (Gl0B  ♦  Gsc)/(B+k)  (6.7-8) 


(a„.  -  v/2tana. . )tana  .s 

B  =  _*i £1 Si-  (6.7-9) 

(na   .    +  v/2tana   .tana., 
si  si  ill 


a..    =   (M..    -  af  -  v2/4)1/2  (6.7-10) 


a  .  =  (M. .  -  a2  -  v2/4)1/2  (6.7-11) 

si     ij    i 


The  modified  eigenvalue  problem  defined  above  is  nonself  adjoint 
and  therefore  its  completeness  is  not  assured.  The  adjoint  problem  will 
have  the  same  form  as  that  obtained  in  Section  3.3. 

Now  if  we  write  the  outer  expansions 


♦  =  I  *k  6k  (6.7-12) 

k=0 


and  collect  the  first  order  terms  in  6,  the  solution  can  be  obtained  in 
the  same  form  as  in  Section  6.3.  For  Se  <  SeoN  we  have 


fj  =  (0,  -lJvDC°,  0,  0)  (6.7-13) 


R  R 

I.  -J  R(3i|)/9r)J  (a.r/R)dr/J  rJ  (a.r/R)dr       (6.7-14) 
1   0         °     l  0   °  l 


h]    =  0  (6.7-15) 

l 
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The  outer  expansion  about  the  planar  state  then  becomes 


"   °°  <Qii'fi>  2 

*  =  ♦  +  S  E   E  — =S — =-  Q.  .J  (a.r/R)  +  0(6  )      (6.7-16) 

i.1  j-1    Mij    «  °  1 


* 

where  Q.  .  is  the  adjoint  eigenf  unction  of  Q.  . .   Similarly  the  outer 
ij  ij 

expansion  about  the  nonplanar  solution  is 


*   1 
*  -  *c  +  e*oNJo(V/R)  +  0(e2)  +  5[Ti^   QN1Jo(V/R)  +  O(e0)]  +  0(*2) 

(6.7-17) 


Once  again  it  can  be  seen  that  both  solutions  will  fail  as 

Se  +  Se  „  and  inner  expansions  are  needed. 
oN  K 


Se(u)  =  Se  „  +  £ub  +  I     5kykb/k!       (6.7-18) 
oN        k=2 


6(m)  =  M°  (6.7-19) 


with  corresponding  expansions  for  the  variables  T.  ,  T  ,  C. ,  C  and  ?. 

X#     5     X*     S 

The  first  order  problem  in  y  then  becomes 


Li*i  =  6'fi  (6-7-20) 
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where  f.  has  been  defined  in  eqn.  (6.7-13).    Here  all  the  boundary 
1 

conditions  will  be  homogeneous  and  the  solvability  condition  gives 


6'  <  qJ1  ,fj>  =  0  (6.7-21) 


Hence  6'  =  0,  C  >  1  and  *  =  A*  .   In  the  next  order 
'  T     o 


L.*2  =  5"f!  (6.7-22) 

iTi       l 


and  at  z  =  0,  B.*2  =  h2  (6.7-23) 


h  =  -2  ASe'  [0,  x,   .  G    ,  0,  o]  +  2A  I..8 
i  L     01   c  11 


The  solvability  condition  now  gives 


So  here  too  we  can  take  6'  '  *  0  and  Se'  *  0 


(6.7-24) 


9     *  9     2 

I.,  =  I  rJ:>(a.r/R)dr/J  rJ  (a.rVR)dr  (6.7-25) 

11    o   °   l       0       °     l 


J(Q*rh2)  -  6"  <  Q^.fJ  >  =  0         (6.7-26) 


Se  -  Se  „  +  561/2  +  0(5)  (6.7-27) 

ON 


*  =  *°  +  A  *  61/2  +  0(6)  (6.7-28) 

T   T     v  o 


where  A  refers  to  A,  and  kot    the  two  solutions  of  the  quadratic  equa- 
v  i  c 

tion  in  A,  eqn.  (6.7-26).   The  positive  solution  A1  will  give  the  inner 
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expansion  matching  the  outer  expansion  eqn.  (6.7-16)  for  Se  <  Se  N  with 
the  outer  expansion  eqn.  (6.7-17)  for  Se  >  SeQN.  The  negative  solution 
A2  will  result  in  the  inner  expansion  matching  eqn.  (6.7-16)  for  Se  > 
SeQN  with  eqn.  (6.7-17)  for  Se  <  SeQN. 

6.8.  Controlling  Imperfections 

In  order  to  calculate  the  imperfection  due  to  heat  loss/gain  the 
form  of  the  functions  f&(z)  and  fg(z)  must  be  known.  A  possible  form 
for  these  are 


fH(z)  =  (Tac(z)  "  V  +  (T1  "  Ta)(z  "  C)/(1  +  °  (6.8-1) 


fsU)  =  (Tsc(z)  "  Ta)  "  T2  "  Ta)(z  "  ?)/(s  "  z)        (6.8-2) 


where  T.  is  the  ambient  temperature  and  T„   and  T   are  T„  and  T  in  the 
A  x.c     sc     x,     s 

planar  state  defined  by  eqns.  (3.2-2)  and  (3.2-3).  This  form  was  chosen 
to  satisfy  the  experimental  conditions  and  to  resemble  Newton's  law  of 
cooling.  Hence  6/k.  and  6/k  will  represent  Biot  numbers  for  the  liquid 
and  solid  sections  of  the  ampoule. 

The  resulting  curves  for  the  two  imperfections  are  shown  in  Figs. 
6-3  and  6-4.  Figure  6-3  corresponds  to  the  case  of  heat  gain  and 
advection  in  the  melt  and  Fig.  6-4  corresponds  to  heat  loss.  In  Chapter 
4  we  concentrated  on  identifying  suitable  growth  conditions  in  order  to 
avoid  morphological  instability.  But  as  we  have  shown  these 
inperfections  undermine  our  earlier  effort  by  causing  the  liquid-solid 
interface  to  be  nonplanar  at  all  times.   All  this  makes  the  situation 
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Se  =  Se0+eSe1 


Se 


Figure  6-3.   Imperfect  bifurcation  diagram  for  heat  loss 
and  advection 


Se  =  Se(J+eSe1 


Figure  6-4.   Imperfect  bifurcation  diagram  for  heat  gain 
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seem  hopeless  for  the  crystal  grower  but  it  should  be  noted  that  even 
though  these  imperfections  cause  nonplanarity  right  from  the  beginning, 
large  deviations  from  planarity  still  occur  only  for  Se  >  Se  „.  Still 
for  those  growers  concerned  even  with  the  small  amount  of  nonplanarity 
at  Se  <  SeoN,  we  suggest  a  possible  way  to  reduce  the  effect  of 
imperfections  by  mutual  cancellation. 

With  regard  to  the  imperfect  insulation  case,  the  important  obser- 
vation to  be  made  is  that  this  is  an  imperfection  that  can  be  controlled 
to  some  extent.  This  can  be  done  either  by  adjusting  the  amount  of 
insulation  or  by  varying  the  ambient  temperature.  If  this  were  the  only 
imperfection  in  the  problem,  its  effect  can  be  minimized  by  good  insula- 
tion and  by  keeping  Ta  between  T1  and  T2.  Coupled  as  it  usually  will  be 
with  the  imperfection  due  to  advection  in  the  melt,  this  control  can  be 
used  to  reduce  their  effects  by  mutual  cancellation.  Total  elimination 
is  of  course  not  possible  but  by  maintaining  the  ambient  temperature 
above  that  of  the  ampoule  it  should  be  possible  to  minimize  the  effect 
of  these  imperfections. 


CHAPTER  7 
NEW  DIRECTIONS 


There  are  several  important  issues  still  waiting  to  be  addressed  in 
morphological  instability  in  crystal  growth.  In  this  Chapter  some  of 
these  which  are  important  and  deserve  attention  will  be  briefly  dis- 
cussed. 

7.1.   Transition  to  Dendritic  Growth 

This  is  a  very  open  question  and  it  will  probably  be  years  before 
even  preliminary  results  can  be  obtained.  Ungar  and  Brown  (1985)  have 
made  an  attempt  in  this  direction  by  modelling  the  formation  of  deep 
cells  but  were  unable  to  proceed  to  the  point  where  the  cells  begin  to 
coalesce.  Another  way  of  looking  at  their  work  is  to  say  that  they 
tried  to  portray  the  transition  to  dendritic  growth  as  a  smooth  one  and 
they  failed  in  this. 

This  brings  up  the  question,  is  the  transition  really  smooth?  If 
it  is  smooth  then  it  can  only  mean  that  the  failure  of  Ungar  and  Brown 
was  only  due  to  shortcomings  of  their  deep  cell  model  and  improving  the 
model  should  rectify  the  problem.  But  if  the  transition  is  similar  to 
transition  to  turbulence  in  fluid  flow  or  transition  to  chaos  in  dynami- 
cal systems,  then  we  have  a  much  more  complex  problem  which  defies  easy 
solutions.  The  experiments  of  Heslot  and  Libchaber  (1984)  indicate  that 
the  coalescing  of  cells  is  catastrophic  rather  than  smooth. 
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A  possible  way  to  determine  whether  the  transition  is  smooth  or  not 
is  to  look  at  secondary  bifurcations  for  the  steady  state  case.  Second- 
ary and  cascading  bifurcations  would  indicate  a  more  complex  transi- 
tion. When  Ungar  and  Brown  (1984a)  initially  did  their  nonlinear  calcu- 
lations with  their  simplified  one-sided  model  they  observed  subcritical 
bifurcations  but  when  they  repeated  their  calculations  for  a  more  accu- 
rate model  (Ungar,  Bennet  and  Brown  (1985))  they  failed  to  see  any. 
Unfortunately  their  calculations  were  performed  for  only  two  sets  of 
experimental  conditions  and  their  results  are  by  no  means  conclusive. 

I  have  already  talked  about  the  multiplicities  in  the  problem  and 
these  could  lead  to  secondary  bifurcations  and  hence  examining  the 
multiplicities  in  the  manner  of  Reiss  (1983)  would  be  a  starting 
point.  If  secondary  bifurcations  are  observed  an  attempt  could  be  made 
to  show  that  dendritic  growth  is  a  result  of  chaos  in  the  system. 

7.2.   Extension  to  Semiconductor  Materials 


The  theory  of  morphological  instability  was  developed  for  binary 
alloys  but  its  biggest  potential  application  is  in  the  growth  of  semi- 
conductor materials  from  the  melt.  But  these  materials  are  different  in 
two  important  respects.  Firstly  they  are  not  simple  binary  mixtures 
like  alloys  but  complex  compounds  in  the  crystal  form  and  secondly  all 
binary  semiconductors  like  gallium  arsenide  have  a  one-to-one  ratio  of 
the  components  and  hence  even  the  melt  is  not  a  dilute  mixture.  The 
effect  of  these  is  that  the  equilibrium  relations  at  the  liquid-solid 
interface  have  to  be  reformulated.  For  example,  the  liquidus  slope  can 
no  longer  be  assumed  to  be  a  straight  line  and  at  least  a  quadratic 
relationship  will  be  required. 
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7.3.  Inclusion  of  Microscopic  Effects 

Until  now  the  researchers  in  morphological  instability  have  been 
fairly  successful  in  studying  the  problem  in  the  presence  of  the  various 
macroscopic  influences  that  beset  it  while  microscopic  effects  have  been 
almost  completely  ignored,  except  for  the  influence  of  grain  boundaries 
(see  Ungar  and  Brown  ( 198Mb)).  Anisotropy,  that  is  the  influence  of  the 
lattice  structure  of  the  crystal  on  morphological  instability,  cannot  be 
ignored  and  this  view  has  been  reinforced  by  the  experiments  of  Heslot 
and  Libchaber  (1984).  But  how  to  include  a  microscopic  effect  on  a 
macroscopic  model? 

This  is  a  fundamental  question  not  just  in  crystal  growth  but  in 
all  fluid  mechanics.  Molecular  dynamics  and  statistical  mechanics  show 
great  promise  in  bridging  the  gap  between  molecular  theory  and  continuum 
mechanics  (see  Bitsanis,  Magda,  Tirrell  and  Davis  (1987)).  Even  though 
this  would  require  tremendous  effort  and  might  seem  like  an  overkill,  it 
is  unavoidable  if  we  hope  to  come  up  with  an  accurate  model  for  the 
influence  of  anisotropy  on  morphological  instability.  More  importantly 
it  will  be  invaluable  in  understanding  the  mechanics  of  solidification 
on  a  molecular  level. 

7.^.  Numerical  Methods 

Almost  all  results  in  this  thesis  have  been  obtained  by  means  of 
asymptotic  expansions.  The  advantages  of  this  have  already  been  dis- 
cussed in  earlier  chapters  but  there  were  numerous  instances  where  the 
effects  of  large  deviations  from  the  basic  problem  are  necessary, 
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requiring  numerical  solutions.   Ideally  both  types  of  analysis  should  be 
done,  one  complementing  the  other. 

Of  course  if  molecular  dynamics  are  employed,  a  numerical  approach 
is  the  only  avenue,  but  even  less  demanding  problems  require  it.  In 
this  work  some  calculations  were  done  for  a  weakly  nonlinear  analysis 
with  hexagonal  and  cylindrical  cells  but  definite  conclusion  about  the 
bifurcation  diagram  could  not  be  reached  due  to  the  absence  of  numerical 
nonlinear  calculations.  Similarly  the  influences  of  large  ampoule  wall 
heat  losses  and  natural  convection  with  rigid  walls  require  numerical 
solutions  not  to  mention  looking  for  secondary  bifurcations  along  the 
nonlinear  curve. 

7.5.   Experiments 

As  mentioned  several  time  in  this  thesis,  there  is  a  strong  need 
for  quantitative  experiments  near  the  onset  points.  Recently  de 
Cheveigne  et  al .  (1985  and  1986)  and  Heslot  and  Libchaber  (1984)  have 
made  some  efforts  in  this  direction  but  a  lot  more  needs  to  be  done  for 
experiments  to  even  catch  up  with  theoretical  work.  Some  of  the  ques- 
tions that  need  to  be  answered  by  experiments  are: 

*  Subcritical  and  supercritical  bifurcations:  As  discussed  in  Section 
4.4  theory  predicts  there  are  subcritical  and  supercritical  bifurca- 
tion regions  for  this  problem  and  experimental  verification  of  these 
would  provide  very  useful  information  for  the  practical  crystal 
grower . 

*  Wavepattern  and  wavenumber  selection:  In  this  work  (as  well  as  in 
the  work  of  others)  wherever  a  choice  had  to  be  made  with  regard  to 
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wavepattern  or  wavenumber,  the  decision  was  based  on  theoretical 
considerations  or  analogy  with  Rayleigh-Benard  convection.  In  some 
instances,  notably  in  wavepattern  selection,  the  issue  was  skirted  by 
considering  all  possibilities.  Experiments  need  to  be  done  to  verify 
the  validity  of  criteria  proposed  or  to  establish  new  criteria  in 
selections. 

*  Reducing  imperfections:  In  Section  6.8  a  way  of  reducing  the  effect 
imperfections  by  playing  off  one  against  the  other  was  proposed. 
Only  experiments  can  tell  if  this  is  a  viable  solution  to  the  problem 
of  imperfections. 

*  Transition  to  dendritic  growth:  Heslot  and  Libchaber  have  shown  that 
oscillatory  solutions  precede  dendritic  growth.  Cell  disappearance 
and  cell  splitting  are  also  reported  precursors.  More  work  is  needed 
to  establish  the  thread  that  leads  from  the  cellular  pattern  near 
onset  to  dendritic  growth. 


APPENDIX  A 
NOMENCLATURE 


Latin 


A  Ratio  defined  in  eqn.  (1.2-24) 

Av  Constant  multiple  defined  by  (6.4-14) 

a  Wavenumber  or  aspect  ratio 

an  Wavenumber  for  n   horizontal  eigenvalue 

a^,  ag  Modified  wavenumbers  defined  in  (3.2-25)  and  (3.2-26) 

a^,  asi  Modified  wavenumbers  defined  in  (5.2-19)  and  (5.2-20) 

a^n  The  wavenumber  for  which  SeQ  is  a  minimum 

B  Boundary  operator 

B  Ratio  defined  in  eqn.  (4.2-23) 

Bi  Biot  number 

Bo  Bond  number 

b  Exponent  defined  by  eqn.  (6.4-1) 

C-j  Boundary  concentration  at  z  =  -I 

C.,C  Solute  concentrations 

c  Exponent  defined  by  eqn.  (6.4-2) 


Cr  Crispation  number 

D 


d_ 
dz 


DN  Constant  used  in  eqn.  (6.4-23) 

Dj,D  Solute  diffusivities 

f  Inhomogeneous  column  vector  in  the  domain 
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f^,  f  Functions  defined  in  (6.2-1)  and  (6.2-2) 

Gj,G  Temperature  gradients  at  the  interface  in  the  planar 

state 
G£c,Gsc  Concentration  gradients  at  the  interface  in  the 

planar  state 

Gc  (V*c  +  DaGsc)/(Dl  +  V 

GT  (kA  ♦  k3Gs)/(kz  +  ka) 

g  Gravitational  vector 

Interface  curvature 

h  Inhomogeneous  column  vector  at  the  boundary 

h  Y/R 

I|  Ratio  of  integrals  defined  by  (6.3-6) 

Ip  Modified  Bessel's  function  of  the  first  kind  of  order 

P 
Ijj  Ratios  of  integrals  for  various  cell  patterns  used  in 

Chapter  1 

J  Bilinear  concomitant 

Jp  Bessel's  functions  of  the  first  kind  of  order  p 

k  Distribution  coefficient 

k„,k  Thermal  conductivities 

%     s 

L  Matrix  differential  operator  for  linear  problem 

L  Latent  heat  number,  L,  D./k.G„A 

h  x,   IT 

L,L,  ,L.2,R  Wavelengths  for  various  horizontal  cell  patterns 

Lh  Latent  heat 

I  Liquid  melt  depth 

I  Total  melt  depth 

M  Eigenvalue  of  augmented  morphological  problem 
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m 


Absolute  value  of  liquidus  slope 
Ma  Marangoni  number 


n 


Unit  normal  vector  at  the  interface,  pointing  into 
the  solid 


p  Constant  defined  in  (4.2-15)  or  modified  pressure 

Pj_.  Ps  Temperature  eigenfunctions  in  augmented  morphological 

problem 

Pe  Peclet  number,  D  /vl 

Pr  Prandtl  number 

Q  Eigenvector  of  augmented  morphological  problem 

q^.  q  Concentration  eigenfunctions  in  augmented 

morphological  problem 

R  Ampoule  radius 

Ra  Rayleigh  number 

s  Solid  depth 

Se  Sekerka  number,  mG  /GT 

^eomin  Tlle  minimura  value  of  Se 

T\,T  ,T  Temperatures 

Ta  Ambient  temperature 

TM  Freezing  temperature  of  pure  solvent 

T-|,T2  Boundary  temperatures  at  z  -  -I   and  z  =  s 

t  Deviation  from  planarity  in  augmented  morphological 

problem 

u  Advective  velocity 

u.,  u  Temperatures  in  inner  expansions 

ur»  uz  Components  of  u 

V  Convective  velocity 

v  Crystal  growth  velocity 
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v.  i  v  Concentrations  in  inner  expansions 

w  Horizontal  factor  of  the  z  component  of  V  or 

deviation  from  planarity  in  inner  expansions 

(kG*c  "  Gsc)/(GZ  "  V 
Y<  Ratio  defined  in  eqn.  (5.2-16) 


Greek 


a,, a  Thermal  diffusivities 

3  ks/k* 

Y  I   /I 

m 

6  Imperfection  parameter 

V„  Surface  Laplacian 

e  Pertubation  variable  about  the  planar  state 

5  Departure  of  interface  from  planarity 

9  Fourier  coefficient  of  temperature 

9.,  6  Transformed  temperatures  defined  by  (6.2-3) 

and  (6.2-4) 

A  Interfacial  tension 

M  61/C 

5  Stretching  parameter  defined  by  eqn.  (6.4-1) 

p0,p  Densities 

a  Eigenvalue  for  linear  stability  problem 

*  Column  vector  (Tn,C„,T  ,C  )  or  (T.,T  ) 

Si   x,   s   s        is 

<J>  n   horizontal  eigenfunction 


Ill 


1  Column  vector  in  inner  expansion 

♦  Adjoint  eigenfunction  of  <j>  for  hexagons 

Q  Inhomogeneous  column  vector  in  eqn.  (6.4-17) 

and  (6.7-24) 


Gothic 


A  Capillary  number,   TJ/LUGJ 

M   hT 


0Q  Modified  concentration  gradient  defined  in  (3.2-28) 

GT  Modified  temperature  gradient  defined  in  (3.2-27) 


Subscripts 

c  Refers  to  the  planar  state 

mn  m  refers  to  order  in  the  perturbation  series  in  e 

n  refers  to  the  nth  horizontal  Fourier  coefficient 
I  Refers  to  liquid  phase 

s  Refers  to  solid  phase 


Superscripts 


Dimensionless  quantity 
Adjoint  vector  or  operator 
Complex  conjugate 
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Variable  in  linear  stability  problem 

Order  of  terms  in  perturbation  series  for  n,  Pe, 

5  or  p 


APPENDIX  B 
PHYSICAL  PROPERTIES 


Properties  of  Alloy  Systems 


Pb-Sn 

Pb-Sb 

C-Austenite 

D  (cm  s   ) 

k.(J  cm   s   K  ) 

k   (J  cm   s   K   ) 
s 

m  (K/wt  %) 

TM  (K) 

-2 
A  (J  cm   ) 

Lh  (J  cm-3) 

k 

-3 
P„  (g  cm  ) 

3 

3.0x10-5 

0.159 

0.297 

2.33 
600 

4.27x10"6 
256 

0.3 
10.66 

2.0x10"5 

0.147 

0.275 
5 
600 

2.30x10"6 
280 

0.4 
10.66 

2.0x10"i* 

0.15 

0.29 
65.3 
1485 

2.04x10"5 
2013 

0.36 

7.4 
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Properties  of  Organic  Systems 


Succi  noni  tr  i  le/ Ace  tone 

CBr1</Br2 

D^    (cm     s      ) 

1  .27x10" 5 

1  .2x10~5 

k^    (J  cm"1    s~1    K~1) 

2.23x10"3 

k      (J  cm       s       K     ) 

2.25x10"3 

1  .1x10"2 

m    (K/wt  %) 

3.02 

2.9 

TM   (K) 

331 

93 

-2 
A    (J  cm     ) 

8.94x10~7 

7.00x10-7 

Lh   (J  cm-3) 

47.8 

3^.45 

k 

0.1 

0.16 

-3 

P_    (g  cm     ) 

s 

1.02 

3.42 
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